CO 

tH 

If) 

tH 

cr 


a: 

LU 

o 


o 

Q 

a: 

LU 


us  Army  Corps 
of  Engineers^ 

Engineer  Research  and 
Development  Center 


INNOVATIVE  SOLUTIONS 
for  a  safer,  better  world 


Operational  Energy  Capabilities  Improvement 

Energy  Performance  Measurement  and 
Simulation  Modeling  of  Tactical  Soft-Wall 
Shelters 

Michael  Deru,  Eric  Bonnema,  Greg  Barker,  Ed  Hancock,  July  2015 

and  Ashok  Kumar 


CD  ^ 
C  O 

LU  O 
E  CD 


O 


O 


CD 

CO  CD 
C  CO 
O  CD 

o  cr 


Approved  for  public  release;  distribution  is  unlimited. 


The  U.S.  Army  Engineer  Research  and  Development  Center  (ERDC)  solves 
the  nation’s  toughest  engineering  and  environmental  challenges.  ERDC  develops 
innovative  solutions  in  civil  and  military  engineering,  geospatial  sciences,  water 
resources,  and  environmental  sciences  for  the  Army,  the  Department  of  Defense, 
civilian  agencies,  and  our  nation’s  public  good.  Find  out  more  at 
www.erdc.usace.army.mil. 


To  search  for  other  technical  reports  published  by  ERDC,  visit  the  ERDC  online  library 
at  http://acwc.sdp.sirsi.net/client/default. 


Operational  Energy  Capabilities 
Improvement 


ERDC/CERLTR-15-13 
July  2015 


Energy  Performance  Measurement  and 
Simulation  Modeling  of  Tactical  Soft-Wall 
Shelters 

Ashok  Kumar 

Construction  Engineering  Research  Laboratory 
U.S.  Army  Engineer  Research  and  Deveiopment  Center 
2902  Newmark  Drive 
Champaign,  iL  61822 

Michael  Deru  and  Eric  Bonnema 

Nationai  Renewabie  Energy  Laboratory 
15013  Denver  West  Parkway 
Goiden,  CO  80401 

Greg  Barker  and  Ed  Hancock 

Mountain  Energy  Partnership 
13900  North  87th  St 
Longmont,  CO  80503 


Final  report 

Approved  for  public  release;  distribution  is  unlimited. 


Prepared  for  Office  of  the  Assistant  Secretary  of  Defense 
3700  Defense  Pentagon 
Washington,  DC  20301-3700 

Under  Project  P2  5D390F,  “Advanced  Energy  Efficient  Shelter  Systems  for 
Contingency  Basing  and  Other  Applications” 


ERDC/CERLTR-15-13 


ii 


Abstract 

Two  soft-wall  tactical  shelters,  Airbeam  Series  32  and  Utilis  TM  60,  were 
set  up  at  an  outdoor  test  facility  in  Champaign,  IL.  Each  was  instrumented 
to  record  energy-loss  measurements  during  a  heating  and  cooling  season, 
and  a  standard  tracer-gas  technique  was  used  to  measure  infiltra¬ 
tion.  These  data,  and  thermal  load  assumptions  for  a  command-operations 
shelter,  were  used  as  inputs  to  develop  and  calibrate  shelter  simulation 
models  using  the  Department  of  Energy’s  EnergyPlus  modeling  platform. 
Simulations  were  run  to  validate  the  calibrated  models,  and  then  a  series 
of  other  simulations  was  run  using  climate  data  from  eight  U.S.  climate 
zones  and  a  locale  in  Kuwait.  These  simulations  included  the  application  of 
energy-conservation  accessories  such  as  shade  flies  and  radiant  heat  barri¬ 
ers  materials. 

Simulation  results  indicate  that  using  either  type  of  radiant  barrier  provid¬ 
ed  the  best  performance  in  cold  climates;  and  the  combination  of  both 
barriers  provided  the  best  performance  in  the  hot  climates.  The  shade  fly 
provided  the  least  savings  in  all  simulated  climates  except  for  Panama 
City,  EL,  where  it  provided  the  highest  savings.  With  limitations  explained 
in  the  report,  the  models  provide  a  useful  technology  for  identifying  energy 
performance  trends  and  making  comparisons  between  the  two  modeled 
shelters. 
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1  Introduction 

1.1  Background 

Current  battlefield  shelter  systems,  including  tents,  are  not  energy  effi¬ 
cient,  and  often  require  powerful,  energy-intensive  heating  and  air  condi¬ 
tioning  units  to  compensate  for  air  leakage.  Expeditionary  base  camps 
consume  excessive  amounts  of  fuel  that  must  continually  be  resupplied, 
which  increases  security  convoy  demand  and  diverts  warfighter  efforts.  A  5 
ton  Environmental  Control  Unit  (ECU)  must  be  transported  to  the  site  of 
each  small  shelter.  There  is  an  urgent  need  to  address  the  energy  demand 
of  forward  operating  bases  (FOBs):  for  every  gallon  of  generator  fuel  con¬ 
sumed  for  sheltering  soldiers,  seven  gallons  are  needed  to  transport  it 
there.  Heating,  ventilating,  and  air  conditioning  (HVAC)  service  comprises 
75%  of  the  electrical  demand,  and  50%  of  that  energy  is  lost  because  soft- 
walled  shelters  are  inefficient  at  retaining  conditioned  air. 

FOBS  are  typically  powered  by  electric  generators  operating  on  diesel  or 
Jet  Propulsion  8  (JP8)  fuel.  These  fuels  are  expensive  and  dangerous  to 
deliver  in  high-risk,  hostile  areas  (Hartranft  2008).  A  large  portion  of  the 
fuel  requirement  is  used  to  condition  the  shelters,  which  provide  little  pro¬ 
tection  against  extreme  environmental  conditions  in  many  areas  of  opera¬ 
tion.  Therefore,  there  is  a  strong  interest  in  finding  materials  and  methods 
that  will  substantially  reduce  energy  consumption  and  improve  soldier 
comfort  inside  shelters. 

The  Engineer  Research  and  Development  Center  (ERDC)  is  a  partner  in  a 
Joint  Service  research  program  supported  by  the  Department  of  Defense 
(DoD)  Operational  Energy  Capital  Improvement  Fund  (OECIF),  which  al¬ 
so  includes  the  Air  Force  Civil  Engineering  Center  (AFCEC)  and  is  led  by 
the  U.S.  Army  Natick  Soldier  Research,  Development,  and  Engineering 
Center  (NSRDEC).  The  program  supports  initial  demonstrations  and  eval¬ 
uations;  development  projects  to  advance  the  state  of  the  art  to  reduce  lo¬ 
gistics  impacts  and  energy  costs;  follow-on  demonstrations  to  optimize 
initial  results  and  validate  cost  savings;  and  transition  efforts  to  field  shel¬ 
ters  that  significantly  reduce  battlefield  fuel  consumption  and  manpower 
requirements. 
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The  demonstration  project  reported  here,  led  by  the  ERDC  Construction 
Engineering  Research  Laboratory  (ERDC-CERL),  investigated  the  energy 
consumption  of  two  soft-sided  shelter  systems  using  an  outdoor  ERDC- 
CERL  FOB  research  facility  and  a  Department  of  Energy  (DOE)  simulation 
modeling  technology. 

1.2  Objectives 

The  objectives  of  this  project  were  to  develop  calibrated  computer  models 
of  two  soft-walled  military  shelters  using  the  EnergyPlus  simulation  mod¬ 
eling  package,  and  then  analyze  those  models  to  investigate  the  energy- 
consumption  impacts  of  adding  insulated  liners,  radiant  barriers,  or  shad¬ 
ing  to  these  shelters  in  different  climates. 

1.3  Approach 

Researchers  applied  data  collection  and  monitoring  methods  developed 
from  research  on  conventional  buildings  to  use  in  calibrating  energy-use 
models  for  two  types  of  soft-sided  military  shelters:  the  HDT  AirBeam 
2032  (HDT  Global,  Solon,  OH)  and  the  Utilis  TM60  (Utilis  USA,  Fort 
Watson  Beach,  FL).  This  report  describes  the  process  of  instrumenting 
and  monitoring  one  shelter  of  each  type  set  up  at  the  ERDC-CERL  For¬ 
ward  Operating  Base  Laboratory  (EFOB-L),  Champaign,  IL,  to  collect  data 
for  input  to  calibrated  baseline  shelter  models  developed  with  the  energy- 
simulation  model  EnergyPlus,  Version  8.0  (DOE  2014a). 

The  baseline  EnergyPlus  models  were  modified  and  analyzed  to  investi¬ 
gate  how  energy  performance  would  be  affected  upon  application  of  vari¬ 
ous  energy-enhancing  accessories  such  as  an  interior  liner,  a  radiant-heat- 
transfer  barrier,  or  a  shade  fly.  These  modifications  were  modeled  in  four 
different  climates.  Also,  the  baseline  energy  model  for  the  Airbeam  shelter 
was  modified  to  investigate  the  performance  of  a  phase-change  material 
(PCM)  liner  in  DOE  climate  zones  1-8. 

The  two  shelters  were  set  up  at  the  ERDC-CERL  Forward  Operating  Base 
Laboratory  (EFOB-L)  facilities  in  Champaign,  IL.  (The  facility  also  in¬ 
cludes  three  hard-sided  shelters  that  are  the  subjects  of  a  separate  study.) 
Field  monitoring  of  the  shelters  was  performed  in  March  and  June  of  2014 
for  cold-  and  hot-weather  testing,  respectively.  The  layout  of  the  shelters  is 
shown  in 
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Figure  i.  The  shading  impact  of  the  B-huts  was  included  in  the  energy 
modeling  of  the  soft-sided  shelters. 


Figure  1.  Site  pian,  June  2014. 


Note:  the  Structurally  Insulated  Panel  Hut  (SIP-Hut)  was  not  in  place  during  the  March  2014  tests. 


The  Utilis  TM6o  is  34  ft  by  19  ft  (10.4  m  by  5.8  m)  by  9  ft  high  (2.7  m), 
with  area  of  650  ft^  (60  m^).  It  is  designed  for  extreme  conditions,  and  can 
be  used  for  various  purposes  such  as  command  and  control  centers,  billet¬ 
ing,  and  medical  (Utilis  2014).  The  shelter,  as  tested,  had  a  Thinsulate™* 
liner  and  a  mesh  shade  fly.  Pictures  of  the  TM60  shelter  with  the  mesh 
shade  fly  set  up  at  the  EFOB-L  are  shown  in  Figure  2. 


*  Thinsulate  is  a  trademark  of  3M,  Maplewood,  MN. 
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Figure  2.  Utilis  TM60  shelter. 


The  HDT  AirBeam  2032  shelter  measures  20  ft  by  32  ft  (6.1  m  by  9.8  m) 
by  11  ft  high  (3.35  m),  with  a  usable  area  of  640  ft^  (59.5  m^)  (HDT  2010). 
The  pressurized  “air  beam”  in  this  shelter  provides  a  wide,  unobstructed 
space  that  can  be  used  for  command  and  control  space,  maintenance  activ¬ 
ities,  or  other  purposes.  The  shelter,  as  tested,  had  a  standard  interior  liner 
and  no  energy-efficiency  features  except  a  vestibule.  Images  of  the 
AirBeam  shelter  with  vestibule  at  ERDC-CERL  are  shown  in  Figure  3. 


Figure  3.  HDT  AirBeam  shelter. 


Each  shelter  was  conditioned  with  a  60,000  Btu  Improved  Environmental 
Control  Unit  (lECU)  from  DRS  Technologies,  Florence,  KY.  This  unit  is 
rated  at  60,000  Btu/h  cooling  capacity  and  30,000  Btu/h  heating  capacity 
at  125  °F  (52  °C)  dry  bulb  (DB)  ambient  and  90  °F  (32  °C)  DB,  and  75  °F 
(24  °C)  wet  bulb  indoor  temperatures.  It  is  rated  at  a  5.8  energy-efficiency 
ratio  at  these  conditions.  The  lECUs  were  only  used  in  the  summer  cooling 
testing.  Electric  resistance  heaters  were  used  during  the  winter  testing. 
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2  Field  Monitoring  and  Data  Collection 

2.1  Overview 

The  general  approach  to  making  measurements  for  the  purpose  of  cali¬ 
brating  a  simulation  model  is  to  simplify  the  system  to  be  calibrated  and 
force  that  system  to  experience  a  large  but  realistic  range  of  conditions 
that  affect  the  behavior  of  the  system.  In  this  case  the  system  is  a  soft- 
walled  shelter,  excluding  the  ground  upon  which  it  sits.  In  order  to  test  the 
shelters  over  a  wide  range  of  operating  conditions  each  shelter  was  tested 
in  a  heating  and  a  cooling  season  (March  and  June,  respectively).  The  tests 
consisted  simply  of  controlling  the  interior  air  temperature  of  each  shelter 
to  a  nearly  constant  value  by  adding  heat  (heating  season)  or  removing 
heat  (cooling  season)  for  a  minimum  of  24  hours,  continuously  measuring 
the  total  amount  of  heat  added  to  or  removed  from  the  shelter.  These  en¬ 
ergy-flow  measurements  were  then  used  along  with  weather,  an  extensive 
temperature  network,  and  other  measurements  to  calibrate  the  energy 
simulation  models  for  the  heating  and  cooling  modes  of  operation. 

Infiltration  and  ground  heat-transfer  energy  flows  are  particularly  difficult 
to  simulate  accurately  in  buildings,  and  even  more  so  in  soft-walled  shel¬ 
ters.  For  calibration  purposes,  these  heat  flows  were  measured  directly  and 
served  as  known  inputs  to  the  calibration  process.  Detailed  measurements 
were  then  used  to  inform  accepted  models  for  infiltration  and  ground  heat 
transfer  for  annual  energy  simulations.  This  approach  allows  for  the  most 
accurate  calibration  feasible  given  the  constraints  of  the  problem  and  the 
need  for  a  reasonably  accurate  simulation  model  If  more  precise  predic¬ 
tions  of  energy  consumption  and  comfort  for  soft-walled  shelters  were  re¬ 
quired,  additional  research  on  infiltration  and  ground  heat  transfer  from 
these  shelters  would  be  needed  for  model  development. 

2.2  Measurements 

Monitoring  plans  were  developed  for  the  heating  season  and  cooling  sea¬ 
son  tests.  The  measurement  points  on  each  shelter  are  specified  in  Table  1 
and  Table  2.  The  layout  of  the  sensors  in  the  shelters  for  the  cooling  season 
test  is  shown  in  Figure  4  and  Figure  5.  The  sensor  placement  for  the  heat¬ 
ing  season  tests  was  the  same  except  for  the  heat  flux  transducers  (HFTs) 
on  the  floor  (which  can  be  seen  in  sec  2. 3. 4.1,  Figure  22  and  Figure  23). 
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Table  1.  Measurement  points  on  each  shelter. 


Measurement 

Instrument 

Quantity 

Outside  fabric  surface  temperature 

Type  T  Thermocouple 

6 

Inside  fabric  surface  temperature 

Type  T  Thermocouple 

6 

Omega  IR  Temperature  Sensor 

4 

Inside  floor  surface  temperature 

Type  T  Thermocouple 

2 

Omega  IR  Temperature  Sensor 

2 

Heat  flux  at  floor 

Hukseflux  HFP03  heat  flux  transducer 

3 

Total  electrical  energy  to  interior  of  shelter 

Continental  Controls  Watt  Node 

1 

Total  electrical  energy  to  ECU  * 

Continental  Controls  Watt  Node 

1 

ECU  supply  air  temperature  &  humidity  * 

VaisalaHMPSO 

1 

ECU  return  air  temperature  &  humidity  * 

VaisalaHMPSO 

1 

^Cooling-season  tests  only 


Table  2.  One-time  measurements  on  each  shelter  (cooling  season). 


Description 

Shelter 

Symbol 

Measurement 

ECU  air  volume  flow  rate 

util  is 

Fv 

42.04  m3/min 

1484.7  cfm 

Air  Beam 

Fv 

40.43  mVmin 

1427.7  cfm 

Equivalent  Leakage 

Area 

Util  is 

ELA 

264  cm2 

41in2 

Air  Beam 

ELA 

466  cm2 

72in2 

Figure  4.  Sensor  piacement  in  the  Utiiis  sheiter. 


-34^25  ft - 


North 


O  =HFT 

O  =  Surface  Temperature 
•  =  Air  Temperature 

O  =  IRTemperature  SenseJ 


■  Sun  Shade 
Notes: 


Sensor#1  = 
SensGr#4  = 


South  End 
North  End 
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Figure  5.  Sensor  placement  in  the  AirBeam  shelter. 


32.0  ft 


North 


T 


1 


O  =HFT 

O  =  Surface  Temperature 
#  =  Air  Temperature 

O  =  IR  Temperature  Sensor 
Notes: 

Sensor  #1  =  South  End 
Sensor #4  =  North  End 


2.2.1  Interior  air  temperature 

Interior  dry-bulb  temperature  was  measured  with  thermocouples  enclosed 
in  mechanically  aspirated  radiation  shields,  shown  in  Figure  6.  The  small 
fan  blows  air  across  the  sensor,  greatly  increasing  the  convection  heat- 
transfer  coefficient  and  reducing  the  effects  of  thermal  radiation  on  the 
sensor.  An  aspirated  radiation  shield  improves  the  accuracy  of  air  temper¬ 
ature  measurement  in  an  environment  where  the  temperatures  of  the  sur¬ 
rounding  surfaces  (wall,  roof,  and  floor)  are  markedly  different  from  the 
air  temperature.  The  interior  dry-bulb  temperature  was  measured  at  two 
points  in  each  shelter  approximately  4  ft  (1.2  m)  above  the  floor,  located 
one-third  of  the  distance  from  each  end  of  the  shelter  and  centered  be¬ 
tween  the  long  sides  of  the  shelter.  The  positions  of  the  sensors  in  the  two 
shelters  are  shown  in  Figure  4  and  Figure  5. 
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Figure  6.  Aspirated  temperature  shieid. 


2.2.2  Surface  temperatures 

Each  shelter  has  an  inner  and  outer  fabric  layer.  These  layers  offer  very 
little  thermal  resistance,  and  the  heat  transfer  between  the  inside  and  out¬ 
side  of  the  shelters  is  dominated  by  convection  and  radiation  at  the  surface 
of  each  layer.  We  placed  thermocouples  in  pairs  (inner  and  outer  surface) 
on  each  of  the  four  sides  of  the  shelters;  and  on  each  of  the  two  sloped  roof 
surfaces  of  the  Utilis  shelter  and  both  sides  of  the  semicircular  cross- 
section  of  the  AirBeam  shelter.  Each  thermocouple  was  attached  to  its  re¬ 
spective  surface  by  adding  a  small  amount  of  thermally-conductive  grease 
to  the  thermocouple  bead  and  taping  the  bead  (plus  6  inches  of  lead  wire) 
using  1  mil  thick  polyester  tape  with  a  silicone  adhesive.  The  thermocou¬ 
ples  on  exterior  surfaces  were  covered  with  a  patch  of  the  shelter  material 
to  provide  the  same  radiative  heat  exchange  as  the  shelter  surface.  These 
surface  temperature  measurements  were  later  compared  to  the  surface 
temperatures  predicted  by  the  simulation  model  as  a  way  of  understanding 
differences  between  the  measured  and  modeled  thermal  performance.  We 
also  placed  two  thermocouples  on  the  top  surface  of  the  floor  (or  insula¬ 
tion  during  the  cooling  season  tests)  in  each  shelter,  in  addition  to  two  in¬ 
frared  temperature  sensors  pointed  at  the  east  and  west  areas  of  the  floor, 
in  order  to  obtain  an  approximation  of  floor  surface  temperature. 
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Figure  7.  IR  temperature  sensors,  with  two  facing 
the  east  and  west  waiis  and  two  facing  the  fioor. 


Figure  8.  Thermocoupie  measuring  outside  surface 
temperature,  covered  with  a  patch  of  the  sheiter  materiai. 


2.2.3  Air  infiltration 

Infiltration— air  exchange  between  the  inside  and  outside  of  the  shelter— is 
a  potentially  significant  portion  of  the  total  heat  transfer  from  a  shelter, 
and  it  is  difficult  to  predict  through  simulation  without  measurements  to 
calibrate  the  models.  The  infiltration  rate  was  measured  directly  during 
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the  testing  period  using  a  tracer  gas  technique  that  introduces  a  small 
amount  of  sulfur  hexafluoride  (SFe)  into  the  shelter  and  measures  its  con¬ 
centration  every  five  minutes.  The  rate  of  decay  in  the  concentration  was 
used  to  calculate  the  rate  of  air  infiltration  in  the  shelter.  This  measured 
infiltration  rate  served  as  a  known  input  to  the  simulation  model  during 
the  calibration  period.  However,  an  improved  model  of  the  infiltration  as  a 
function  of  the  driving  functions  (wind  speed,  wind  direction,  and  in¬ 
side/outside  temperature  difference)  is  needed  for  better  predictions  of 
the  infiltration  during  the  annual  simulations. 

For  simple,  single-zone  buildings,  EnergyPlus  estimates  the  infiltration 
with  an  approach  based  on  the  Sherman-Grimsrud  model  (DOE  2014b),  as 
shown  in  Equation  1.  Time-dependent  inputs  to  the  model  include  wind 
speed  and  inside/outside  temperature  difference,  while  constants  affecting 
the  prediction  include  the  equivalent  leakage  area  (ELA)  and  coefficients 
for  the  temperature  and  wind-speed  terms.  The  schedule  multiplier  can  be 
used  to  model  different  modes  of  operation  that  affect  the  infiltration  level. 

Infiltration  =  ■  AT  +  -V^)  (1) 


where 

Fsch 

ELA 

AT 

Fu. 

Cs 


scheduled  fractional  multiplier 
equivalent  leakage  area  at  4  Pa  [cm^] 
average  inside-outside  temperature  [°C] 
average  wind  speed  [m/s] 

stack  or  temperature  coefficient  [(L/s)2/(cm4-K)] 
wind  pressure  coefficient  [(L/s)2/(cm4-(m/s)2)] 


The  ELA  is  most  commonly  measured  by  performing  a  blower-door  test  on 
the  building,  where  the  blower  is  used  to  either  pressurize  or  depressurize 
the  building  over  a  range  of  flow  rates,  resulting  in  a  range  of  in¬ 
side/outside  pressure  differences.  The  flow  versus  pressure  difference 
curve  is  best-fit  to  the  form 


F^  =  C-  AP” 


(2) 


where 


Fv  =  volume  flow  rate 
AP  =  inside-outside  pressure  difference 
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C,n  =  constants  derived  from  a  regression  of  the  data 

The  blower-door  fan  was  installed  in  the  rigid  doorway  of  the  Utilis  shelter 
as  shown  in  Figure  lo,  and  pressurization  tests  were  performed  over  a 
range  of  15  -  35  Pa.  Blower-door  tests  in  buildings  are  performed  up  to  50 
Pa,  but  we  were  concerned  that  the  shelter  might  be  damaged  under  too 
much  pressure  (although  a  50  Pa  pressure  differential  can  certainly  be  in¬ 
duced  by  a  strong  wind).  The  resulting  flow-rate  curve  is  shown  in  Figure 
9,  and  the  ELA  was  determined  to  be  41  in^  (264  cm^). 


Figure  9.  Utilis  shelter  blower-door  test  results. 


Figure  10.  Blower  door  installed  on  the  Utilis  shelter. 
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The  AirBeam  shelter  does  not  have  a  rigid  doorway  to  mount  the  blower- 
door  fan,  so  we  created  a  frame  for  the  fan  out  of  polystyrene  insulating 
board  and  taped  it  into  the  south  doorway  as  shown  in  Figure  ii.  On  the 
AirBeam  shelter  we  completed  both  pressurization  and  depressurization 
tests  from  lo  -  30  Pa.  The  results  of  the  two  tests  are  shown  in  Figure  12, 
and  the  ELA  from  the  pressurization  tests  was  determined  to  be  72  in^ 
(466  cm2). 


Figure  11.  Blower-door  fan  installed  in  the  AirBeam  shelter. 


Figure  12.  AirBeam  shelter  blower-door  test  results. 
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The  equation  fit  for  the  depressurization  tests  in  Figure  12  is  not  as  good  as 
for  the  pressurization  tests,  especially  the  three  points  at  the  highest  pres¬ 
sures.  We  believe  that  this  occurred  because  the  conditions  were  not  at  a 
steady  state  while  the  measurements  were  being  taken;  in  fact,  we  ob¬ 
served  that  the  floor  of  the  shelter  was  gradually  rising  while  the  shelter 
was  under  negative  pressure  (see  Figure  13).  When  inflating,  the  shelter 
comes  to  a  steady  state  fairly  quickly  (on  the  order  of  30  seconds);  howev¬ 
er,  we  never  saw  it  come  to  a  steady  state  while  deflating  because  the  vol¬ 
ume  of  the  shelter  was  still  changing.  For  comparison,  a  blower-door  test 
on  a  rigid  structure  quickly  comes  to  a  steady  state  condition. 


Figure  13.  The  floor  of  the  AirBeam  shelter 
gradually  rising  during  the  depressurization  test. 


2.2.4  Heating  and  cooling  energy 

The  heating-  and  cooling-season  tests  were  performed  by  inducing  a  large 
heat  flow  on  the  shelters  that  could  be  measured  accurately  over  the  peri¬ 
od  of  the  test.  For  each  test,  the  shelters  were  conditioned  to  constant 
temperatures  to  create  a  large  temperature  difference  across  the  envelope 
(>  30  °F).  If  the  shelters  were  simply  allowed  to  “float”  with  little  or  no  ar¬ 
tificially  induced  heat  flow,  it  would  be  difficult  to  adjust  any  assumptions 
in  the  model  to  match  the  measurements.  For  example,  with  very  little 
temperature  difference  between  the  inside  and  outside  of  the  shelter,  the 
simulated  thermal  conductivity  of  the  walls  could  be  changed  by  a  large 
factor  without  affecting  the  simulated  interior  temperature. 

2.2.4.1  Heating  season 

The  shelters  were  heated  using  electric  convection  heaters  to  a  tempera¬ 
ture  of  64.4  °F  (18  °C)  for  the  duration  of  the  test.  This  approach  provides 
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a  well  known  energy  input  and  a  strong  driver  to  force  a  large  heat- 
transfer  rate  from  the  inside  to  the  outside  of  the  shelter.  Small  fans  were 
installed  in  each  shelter  to  attempt  to  keep  the  air  inside  the  shelter  well 
mixed  in  order  to  simplify  the  assumptions  about  stratification  in  the  sim¬ 
ulation  model.  Total  electrical  power  entering  the  shelter  (heaters,  fans, 
lights,  data  acquisition  equipment)  was  measured  for  each  shelter.  All  of 
this  electrical  energy  eventually  becomes  sensible  internal  gain  to  the  shel¬ 
ter  and  was  included  in  calibrating  the  model. 

2. 2. 4. 2  Cooling  season 

The  shelters  were  cooled  using  the  standard  heat  pump  units  used  by  the 
Army,  known  as  ECUs  (Environmental  Control  Units),  to  a  controlled 
temperature  of  approximately  64.4  °F  (18  °C)  in  the  Utilis  shelter  and  68 
°F  (20  °C)  in  the  AirBeam  shelter  through  the  duration  of  the  test.  The  set 
points  were  arbitrary;  in  fact,  the  thermostats  on  the  ECUs  have  no  values 
written  on  them,  just  an  analog  hot-to-cold  dial.  The  goal  was  simply  to 
induce  a  large  heat-transfer  rate  from  the  (hot)  outside  to  the  (cool)  inside 
of  the  shelter. 

Because  the  ECUs  deliver  air  to  the  shelters  through  a  series  of  air  jets  well 
distributed  along  the  length  of  the  shelter,  it  was  expected  that  the  air  in 
the  shelters  would  be  isothermal.  Total  electrical  power  entering  the  shel¬ 
ter  (lights,  data  acquisition  equipment,  etc.)  was  measured  for  each  shelter 
and  added  to  the  calibration  model  as  sensible  internal  gain  to  the  shelter. 

The  normal  operation  mode  of  an  ECU  is  with  constant  evaporator  fan  op¬ 
eration,  and  compressor  and  condenser  fan  engaging  as  cooling  is  re¬ 
quired.  The  air  volume  flow  rate  of  each  ECU  was  measured  using  a 
“TrueFlow  Plate”  manufactured  by  The  Energy  Conservatory,  Minneap¬ 
olis,  MN  (see  Figure  14).  This  is  a  device  with  a  number  of  calibrated  ori¬ 
fices  and  an  array  of  pitot  tubes  used  to  measure  the  air-flow  rate  in 
heating  and  cooling  systems,  and  has  a  published  accuracy  of  +/-  7%.  We 
assumed  that  the  ECU  air  flow  rate  was  maintained  relatively  constant 
throughout  the  testing  period. 
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Figure  14.  Temporary  installation  of  the  TrueFlow 
Plate  for  a  one-time  measurement  of  ECU  air  flow. 


T7 
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Temperature  and  humidity  sensors  were  placed  in  the  air  streams  of  both 
the  supply  and  return  ducts  in  each  shelter  as  shown  in  Figure  15.  Using 
the  one-time  measurement  of  air  flow  rate  and  the  measured  temperature 
and  humidity  conditions  at  the  supply  and  return  ducts,  total,  sensible, 
and  latent  heat  removed  by  the  ECU  were  calculated  in  real  time  in  the  da¬ 
ta  acquisition  system. 


Figure  15.  Measurement  point 
for  return  air  temperature  and  humidity. 
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2.2.5  Floor  surface  heat  flux 

Heat-transfer  rate  to  the  ground  was  measured  using  three  heat  flux 
transducers  in  each  shelter.  Application  of  the  transducers  was  slightly  dif¬ 
ferent  between  heating-  and  cooling-season  tests.  The  measured  heat  flux¬ 
es  through  the  floor  were  used  as  inputs  to  the  calibration  models  with  no 
attempt  to  model  the  behavior  of  the  ground  itself. 

2.2. 5.1  Heating  season 

Heat-transfer  rate  to  the  ground  was  measured  using  three  heat  flux 
transducers  in  each  shelter.  The  heat  flux  transducers  were  placed  on  the 
floor  with  heat-conductive  grease  to  ensure  good  thermal  contact  with  the 
shelter  floor  fabric.  Sensor  i  was  placed  in  the  center  of  the  floor;  sensor  3 
was  placed  approximately  6  in.  (15  cm)  from  the  outer  edge  of  the  floor; 
and  sensor  2  was  placed  at  a  point  equidistant  between  the  first  two  (see 
section  2.3.4.1,  Figure  22  and  Figure  23). 

2.2.5.2  Cooling  season 

The  heat  flux  transducers  were  adhered  to  the  polystyrene  floor  insulation 
with  heat-conductive  grease  to  ensure  good  thermal  contact  with  the  poly¬ 
styrene  sheet.  The  sensors  were  placed  at  the  center  of  the  floor  on  the  top 
surface  of  the  polystyrene  (sensor  2),  in  the  center  of  the  floor  on  the  bot¬ 
tom  surface  of  the  polystyrene  (sensor  1),  and  at  a  point  about  12  in.  (30 
cm)  from  the  eastern  edge  on  the  top  surface  of  the  polystyrene  (sensor  3), 
as  shown  previously  in  Figure  4  and  Figure  5.  Sensors  1  and  2  were  placed 
on  the  top  and  bottom  of  the  insulating  sheet  in  the  same  central  location 
to  confirm  that  the  two  sensors  read  the  same  heat  flux  when  averaged 
over  time;  it  was  expected  that  the  sensor  on  the  top  would  react  quickly  to 
changes  in  convective  and  radiative  heat  transfer,  while  the  one  on  the  bot¬ 
tom  would  be  a  damped  signal  representing  the  signal  from  the  top  sensor 
averaged  over  time. 

2.2.6  Weather 

Outdoor  weather  conditions  of  dry-bulb  temperature,  relative  humidity, 
wind  speed,  wind  direction,  and  global  horizontal  solar  radiation  were 
measured  during  the  tests  using  field-grade  instruments  mounted  on  a 
temporary  tripod  near  the  shelters.  These  measurements  were  used  as 
guides  during  the  test,  but  for  comparison  of  model  to  field  data,  meas¬ 
urements  collected  onsite  concurrently  by  Oak  Ridge  National  Laborato- 
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ries  (ORNL)  were  used.  An  ORNL  weather  station  set  up  to  support  a  sep¬ 
arate  EFOB-L  project  included  research-grade  instruments  for  measuring 
air  temperature  and  humidity,  wind  speed  and  direction,  beam  and  diffuse 
components  of  solar  radiation,  and  infrared  sky  radiation.  The  direct 
measurement  of  sky  radiation  was  used  as  a  direct  input  to  the  simulation 
model,  bypassing  the  correlation  typically  used  in  building  simulations  for 
predicting  sky  radiation. 

2.3  Observations  for  model  calibration 

2.3.1  Infiltration 

The  infiltration  rate  calculated  from  the  decay  rate  of  SF6  in  each  shelter 
varied  from  o.i  to  0.6  air  changes  per  hour  (ACH).  As  we  observed  during 
our  tests  in  March,  the  infiltration  rate  appeared  to  be  not  only  highly  cor¬ 
related  with  wind  speed  and  inside/outside  temperature  difference  (as  ex¬ 
pected),  but  with  wind  direction.  When  the  wind  was  coming  from  the 
northwest,  the  AirBeam  shelter  showed  much  higher  infiltration  rates  than 
the  Utilis.  When  the  wind  was  in  the  opposite  direction,  the  opposite  trend 
was  true.  Two  possible  explanations  for  this  are  that  (i)  the  Utilis  has  its 
door  facing  south  while  the  AirBeam  has  its  door  facing  north,  and  (2)  the 
Utilis  is  sheltered  by  a  building  to  its  northeast  and  the  AirBeam  is  shel¬ 
tered  by  a  building  to  its  south. 

This  behavior  illustrates  the  need  to  use  measured  infiltration  rates  as  an 
input  to  the  simulation  model  in  order  to  calibrate  the  model  of  the  shel¬ 
ters.  Trying  to  predict  the  infiltration  rates  under  the  exact  conditions  at 
FRDC-CFRL  during  the  tests  would  likely  prove  unsuccessful  because  the 
typical  approach  in  building  simulation,  using  the  Sherman-Grimsrud 
model  (DOE  2014b)  and  blower-door  results,  does  not  account  for  small 
wind  disturbances  that  appear  to  play  a  major  role  in  the  infiltration  rates 
of  the  shelters.  Results  from  the  blower-door  test  were  used  to  generate 
inputs  to  models  for  the  purpose  of  long-term  prediction  of  shelter  infiltra¬ 
tion  rates. 

Because  the  concentration  of  the  tracer  gas  introduced  in  each  shelter 
gradually  diminishes  over  time  as  the  indoor  air  is  mixed  with  outdoor  air, 
the  concentration  eventually  becomes  close  to  the  outdoor  background 
level,  at  which  point  the  infiltration  rate  cannot  be  inferred  from  the  con¬ 
centration  measurement.  When  the  infiltration  rate  was  high  (near  1  ACH) 
the  total  time  between  the  initial  injection  of  the  gas  and  the  point  at 
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which  the  concentration  was  too  low  to  measure  was  on  the  order  of  5 
hours.  Because  the  research  team  did  not  have  access  to  the  site  between 
1700  and  0500  hours  the  following  day,  there  were  significant  periods 
during  which  the  concentrations  of  tracer  gas  were  too  low  and  infiltration 
measurement  rates  could  not  be  collected.  To  fill  in  the  gaps  in  the  data, 
we  created,  correlations  to  fit  the  measured  infiltration  rates  to  the  meas¬ 
ured  wind  speed  and  direction  (as  noted  earlier,  the  wind  direction  ap¬ 
peared  to  have  a  large  impact  on  the  infiltration  rate  of  each  shelter),  and 
inside/outside  temperature  difference.  The  three  independent  variables  in 
the  fit  were  defined  as  wind  speed  (m/s),  the  cosine  of  the  wind  direction 
in  degrees,  and  inside/outside  temperature  difference.  The  cosine  was 
chosen  to  allow  the  wind-direction  variable  to  have  a  continuous  function; 
in  contrast,  the  wind  direction  in  degrees  is  a  discontinuous  function 
where  it  transitions  from  359.999  degrees  (almost  exactly  north)  to  0.000 
degrees  (exactly  north).  A  sine  function  would  have  served  equally  well. 
Equation  3  is  the  best-fit  equation  for  the  Utilis  measured  infiltration  data 
and  Equation  4  was  the  best  for  the  AirBeam,  both  during  the  heating- 
season  tests.  The  forms  of  the  best-fit  equations  have  no  particular  physi¬ 
cal  significance.  The  uncertainties  in  the  fits  are  illustrated  by  graphing  the 
predicted  infiltration  rates  versus  the  measured  infiltration  rates  in  Figure 
16  and  Figure  17  for  the  two  shelters. 

ACH  =  a  +  bx^  +  CX2  -I-  dxl  +  exf  -I-  fxiX2  +  gxf  +  /ixf  -I-  -I-  jxlx2  (3) 

ACH  =  ax^c^^  (4) 


where 

Xi  =  wind  speed  [m/s] 

X2  =  cos  (wind  direction  [deg]) 
X3  =  AT  [oC] 


ACH  (Predicted)  |  ACH  (Predicted) 
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Figure  16.  Goodness  of  fit  of  air-change  rates  as  a  function  of 
wind  speed,  wind  direction,  and  inside/outside  temperature 
difference  for  the  Utiiis  sheiter  in  the  heating  season. 
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Figure  17.  Goodness  of  fit  of  air-change  rates  as  a  function  of 
wind  speed,  wind  direction,  and  inside/outside  temperature 
difference  for  the  AirBeam  sheiter  in  the  heating  season. 
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2.3.2  Floor  heat  flux 

From  experience  with  buildings,  we  expect  the  heat  flow  through  the  floor 
to  be  largest  near  the  edge  since  that  represents  a  shorter  heat-flow  path  to 
the  outdoors  and  the  ground  under  the  center  of  the  building  is  closer  to 
the  building  temperature.  We  also  normally  expect  the  heat-transfer  rate 
at  any  point  on  the  floor  to  be  fairly  steady  if  the  indoor  air  temperature  is 
steady  because  buildings  have  thermally  massive  layer  in  contact  with  the 
ground  that  dampens  variations  in  the  heat  flows.  With  the  shelters  we 
found  neither  of  these  expectations  to  be  correct.  Likely  explanations  in¬ 
clude  the  thin  floor  liner  in  contact  with  the  ground,  the  ground  under  the 
shelter  was  close  to  the  undisturbed  ground  surface  temperature,  and  that 
the  heat  transfer  to  the  ground  in  the  shelters  was  very  dependent  on  con¬ 
vective  heat  transfer  with  the  air  and  radiant  exchange  with  the  inside 
shelter  surfaces,  whose  temperatures  changed  dramatically  depending  on 
the  outdoor  conditions.  In  a  typical  insulated  building  the  inside  surface 
temperatures  do  not  change  much  due  to  a  change  in  the  outdoor  condi¬ 
tions. 

2.3.2.1  Heating  season 

Because  these  shelters  were  unconditioned  until  we  began  heating  them 
for  our  tests,  the  ground  started  out  cold  and  was  expected  to  be  a  large 
heat  sink.  We  placed  3  heat  flux  transducers  on  the  floor  of  each  tent,  one 
in  the  center,  one  near  the  edge,  and  one  between  the  first  two. 

The  hourly-averaged  heat  flux  measurements  for  each  shelter  are  shown  in 
Figure  18.  Note  that  the  trend  among  the  three  measurements  is  different 
between  the  two  shelters.  The  highest  heat-transfer  rate  in  the  Utilis  is 
near  the  edge  as  measured  by  heat  flux  transducer  3  (HFT  3),  but  the 
highest  rate  in  the  AirBeam  is  at  the  center  (HFT  1).  The  difference  is  likely 
due  to  the  radically  different  radiant  environment  created  by  the  insulated 
surfaces  of  the  Utilis  compared  with  the  conditions  created  by  the  uninsu¬ 
lated  surfaces  of  the  AirBeam. 
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2.3. 2.2  Cooling  season 

For  the  summer  cooling  season  monitoring,  the  Utilis  shelter  was  condi¬ 
tioned  (cooled)  until  a  few  days  before  testing  began,  while  the  AirBeam 
shelter  had  been  unconditioned  for  many  months.  The  ground  tempera¬ 
ture  was  an  unknown  initial  condition  for  both  shelters;  however,  it  was 
expected  to  be  cooler  for  the  Utilis  shelter  because  of  the  conditioning,  and 
was  expected  to  be  warmer  for  the  AirBeam  shelter.  To  reduce  the  magni¬ 
tude  of  the  heat  flow  through  the  floor  we  placed  i  in.  thick  extruded  poly¬ 
styrene  insulation  on  the  floors  of  both  shelters.  We  placed  three  HFTs  on 
the  floor  of  each  shelter,  as  discussed  in  section  2.2.5. 

The  15  minute  averaged  floor  heat  flux  measurements  for  each  shelter  are 
shown  in  Figure  19.  The  trends  among  the  three  measurements  are  similar 
between  the  two  shelters;  the  top-surface  heat  fluxes  are  similar  in  magni¬ 
tude  regardless  of  whether  they  are  close  to  the  edge  or  in  the  center,  and 
their  values  change  rapidly  on  a  15  minute  basis  compared  to  the  heat 
fluxes  measured  underneath  the  insulation.  As  expected,  the  rapidly- 
changing  values  measured  at  the  top  surface  of  the  insulation  average  out 
to  the  slowly-changing  values  measured  at  the  bottom  surface  of  the  insu¬ 
lation.  A  trend  that  is  different  between  the  two  shelters  is  that  the  overall 
direction  of  the  heat  flux  in  the  Utilis  shelter  was  from  the  ground  to  the 
inside  of  the  shelter,  while  for  the  AirBeam  shelter  the  heat  flux  tended  to 
be  from  the  shelter  to  the  ground  during  the  middle  of  the  day.  This  is  like¬ 
ly  because  the  AirBeam  had  no  insulating  liner  so  that  when  the  roof  and 
wall  surfaces  of  the  shelter  became  hot  due  to  solar  radiation,  the  heat  was 
radiated  to  the  cooler  ground  surface.  The  insulating  liner  in  the  Utilis 
shelter  greatly  attenuated  this  effect. 


Infrared  images  were  taken  during  the  heating  season  testing  to  provide  a 
qualitative  understanding  of  the  heat-transfer  paths  and  relative  magni¬ 
tudes  in  the  shelters.  The  images  in  Figure  20  show  the  high  heat  flow 
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around  the  door,  at  the  edge  of  the  floor-wall  intersection,  and  through  the 
points  where  the  Thinsulate  liner  is  tacked  together  in  the  Utilis  shelter. 
The  image  of  the  floor  also  shows  the  increased  heat  flow  where  the  floor 
liner  is  in  contact  with  the  ground.  The  images  of  the  AirBeam  shelter  in 
Figure  21  show  the  high  flows  through  the  door  and  the  window.  The  win¬ 
dow  is  closed  in  the  image,  but  the  heat  flow  is  higher  here  because  there  is 
no  air  gap  between  the  inner  liner  and  the  outer  shell. 


Figure  20.  Infrared  images  of  the  Utilis  shelter  door,  walls,  and  floor. 


Figure  21.  Infrared  images  of  the  AirBeam  shelter  door,  walls,  and  window. 


2.3.4  Total  heat-flow  calculations 

2. 3.4.1  Floor  heat-flow  calculations  (Qt) 

A  weighted  average  of  the  HFT  measurements  was  calculated  to  represent 
the  total  heat  flux  through  the  floor  for  both  the  heating-  and  cooling- 
season  tests.  Weights  were  based  on  the  floor  areas  shown  in  Figure  22 
and  Figure  23. 
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Figure  22.  AirBeam  floor  area  weighting  for  heat  flux  measurements. 
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Figure  23.  Utilis  floor  area  weighting  for  heat  flux  measurements. 
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Because  the  floors  were  not  insulated  during  the  heating-season  tests,  the 
total  ground  heat  flux  was  a  significant  portion  of  the  total  heating  load  of 
the  shelter.  It  was  expected  that  the  measured  heat  fluxes  would  be  an  ex¬ 
aggerated  value  since  the  HFTs  were  intentionally  held  down  against  the 
floor  with  no  air  gap  beneath  them.  From  a  visual  inspection,  it  was  esti¬ 
mated  that  approximately  half  the  floor  area  had  small  gaps  or  bubbles  of 
air  between  the  floor  fabric  and  the  ground.  As  a  best-guess  approach  to 
calculating  the  total  heat  flux  through  the  floor  for  the  heating-season 
tests,  we  assumed  that  50%  of  the  floor  fabric  was  touching  the  ground 
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while  the  other  half  had  a  i  cm  (0.4  inch)  air  gap  beneath  it.  It  was  further 
assumed  that  the  temperature  difference  observed  between  the  warmer 
(with  air  gaps)  and  cooler  (without  air  gaps)  sections  of  the  floor  was  ap¬ 
proximately  3.6  °F  (2  °C).  This  value  was  inferred  by  observing  the  tem¬ 
perature  differences  shown  in  infrared  photographs  of  the  floor  in  Figure 
20. 


Total  heat  flux  through  the  floor  was  calculated  using  Equation  5: 


<?/-(!“  fair)  (^iQl  +  A2Q2+  A2Q3)/ (Ai  +  A2+  A3)  -I-  fairiffair/ d)AT )  (5) 


where 


Of 

fair 

Ai^2^3 

QhQ2,Q3 

kair 

d 

AT 


total  heat  flux  through  the  floor  [W/m^] 
fraction  of  floor  with  air  gaps 
weighting  areas  [m^] 
heat  flux  measurements  [W/m^] 
thermal  conductivity  of  air  [W/m-K] 
thickness  of  air  gap  [m] 

air  gap  temperature  difference  between  fabric  and  ground  [K] 


In  calculating  Qf  for  the  cooling  season  data,  fair  was  assumed  equal  to  zero 
since  the  floors  were  covered  in  rigid  insulation  and  any  air  gaps  below  the 
floor  would  have  an  insignificant  effect  on  the  uniformity  of  heat  flux 
across  the  floor. 


2.3.4.2  Infiltration  thermal  load  (Qinf) 

Heat  lost  due  to  infiltration  was  calculated  according  to  Equation  6: 

Qinf  =  ACH  ■  V  ■  Cp  ■  p  ■  (Tin  -  Tout)  (6) 

where 

Qinf  =  rate  of  energy  lost  due  to  infiltration  [W] 

ACH  =  rate  of  infiltration,  air  changes/hour 

V  =  interior  volume  of  shelter  (cross-sectional  area  times  length) 
[m3] 

Cp  =  specific  heat  of  outdoor  air  [J/kg-K] 
p  =  density  of  outdoor  air  [kg/m3] 

Tin  =  average  interior  dry  bulb  temperature  [°C] 
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Tout  =  outdoor  dry  bulb  temperature  [°C] 

2.3.4.3  Mechanical  heating  or  cooling  total  heat  flow  (Qtot) 

For  the  heating-season  tests,  the  total  heat  flow  due  to  mechanical  heating 
was  the  total  electrical  power  to  all  devices  within  the  shelter.  For  the  cool¬ 
ing-season  tests,  the  total  heat  flow  into  the  shelter  was  calculated  as  the 
energy  added  by  all  electrical  devices  within  the  shelter  minus  the  total 
sensible  energy  removed  by  the  ECU,  as  shown  in  Equation  7: 


Qtot  ~  Qelec  ^Cp{Tret  ^sup) 


(7) 


where 

Qtot  =  total  sensible  energy  added  to  the  shelter  via  mechanical 
means  [W] 

Qelec  =  total  electrical  energy  added  to  interior  of  shelter  [W] 
m  =  mass  flow  rate  of  air  through  ECU  [kg/s] 

Cp  =  specific  heat  of  air  [J/kg-K] 

Tret  =  ECU  return  air  temperature  [°C] 

Tsup  =  ECU  supply  air  temperature  [°C] 

Equation  7  is  valid  for  both  heating  and  cooling  seasons,  where  the  ECU 
air  mass  flow  rate  was  equal  to  zero  during  the  heating-season  tests.  Eor 
the  cooling  season  Qtot  was  generally  negative,  indicating  more  energy  re¬ 
moved  from  the  shelter  than  added  to  it  by  electrical  devices  within  the 
shelter. 

2.3.4.4  Wall  and  roof  total  heat  flow  (Qwr) 

The  total  heat  flow  through  the  walls  and  roof  of  each  shelter  was  calculat¬ 
ed  using  Equation  8: 


Qwr  ~  Qtot  Qf^f  Qinf 


(8) 


where 


Af  =  floor  area  [m^] 


ERDC/CERL  TR-15-13 


28 


The  calculated  components  of  heat  flow  for  each  shelter  for  the  heating- 
and  cooling-season  tests  are  shown  in  Figure  24  through  Figure  27.  A  few 
trends  are  noteworthy: 

•  The  cooling  load  due  to  ground  heat  flow  was  almost  negligible  during 
the  cooling-season  tests,  whereas  it  was  a  significant  portion  of  the  to¬ 
tal  heating  load  for  the  heating-season  tests.  This  was  largely  because 
we  had  installed  insulation  on  the  floors  of  both  shelters  for  the  cool¬ 
ing-season  tests  to  minimize  the  ground  heat  flow  and  thus  minimize 
the  importance  of  uncertainty  in  its  measurement. 

•  The  magnitudes  of  the  heat  flows  through  the  walls  and  roof  of  each 
shelter  were  much  smaller  for  the  Utilis  shelter  than  for  the  AirBeam 
shelter.  Some  of  this  effect  was  due  to  the  larger  size  of  the  AirBeam 
shelter,  but  the  most  significant  cause  was  most  likely  that  the  Utilis 
shelter  had  an  insulated  liner  installed. 

•  The  magnitudes  of  the  heat  transfer  due  to  infiltration  were  much 
smaller  during  the  cooling-season  tests  for  both  shelters.  This  was  due 
to  a  combination  of  lower  wind  speeds  and  lower  indoor/outdoor  tem¬ 
perature  differences  during  the  cooling-season  tests. 


Figure  24.  AirBeam  winter  heat  flows. 
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Figure  25.  Utilis  winter  heat  flows. 
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Figure  27.  Utilis  summer  heat  flows. 
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3  Model  Development  and  Calibration 

3.1  Method 

EnergyPlus  Version  8.i  (DOE  2014a)  was  used  to  develop  the  energy  mod¬ 
els  of  the  AirBeam  and  Utilis  shelters.  These  models  were  developed  based 
on  engineering  drawings,  envelope  material  specifications,  material  prop¬ 
erty  measurements,  ECU  specification  sheets,  field  measurements,  and 
observations. 

3.1.1  Geometry 

The  geometry  of  the  AirBeam  and  Utilis  shelters  are  shown  in  Figure  28 
and  Figure  29,  respectively.  The  air  gap  formed  by  the  inner  and  outer  lin¬ 
ers  was  represented  by  individual  zones  surrounding  the  inner  occupied 
zone.  Mixing  objects  between  air-gap  zones  were  used  to  model  the  air 
flow  in  the  air  gap.  The  air  flow  in  this  gap  was  set  to  a  constant  value  of 
about  1 CFM  (0.00049  m3/s). 

Figure  28.  AirBeam  sheiter  modei. 
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Figure  29.  Utilis  shelter  model. 


3.1.2  Material  properties 

The  surface  properties  of  the  inner  shell,  outer  liner,  and  shade  fly  were 
based  on  laboratory  measurements  performed  at  NREL.  Thermal  re¬ 
sistance  properties  came  from  manufacturer  data  or  estimates  based  on 
engineering  tables  for  similar  materials.  The  properties  are  given  in  Table 
3  and  Table  4. 


Table  3.  AirBeam  shelter  material  properties  for  model  calibration. 


Property 

Units 

Exterior  Shell 

Inner  Liner 

Floor 

Outer 

Surface 

Inner 

Surface 

Outer 

Surface 

Inner 

Surface 

Inner 

Surface 

Roughness 

- 

Medium 

Smooth 

Medium 

Smooth 

Medium 

Smooth 

Medium 

Smooth 

Medium 

Smooth 

Thermal  absorptance 

Fraction 

0.897 

0.90 

0.889 

0.899 

0.72  V  0.742b 

Solar  absorptance 

Fraction 

0.571 

0.9 

0.13 

0.13 

0.2 

Visible  absorptance 

0.6 

0.9 

0.13 

0.13 

0.2 

Thermal  resistance 

m2.K/W 

0.0088 

0.0088 

0.5V  0.0176b 

a  Surface  property  of  polystyrene  for  cooling  season  tests 
b  Surface  property  of  standard  floor  material  for  heating  season  tests 
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Table  4.  Utilis  shelter  material  properties  for  modei  caiibration. 


Property 

Units 

Exterior  Shell 

Inner  Liner 

Floor 

Outer 

Surface 

Inner  Surface 

Outer 

Surface 

Inner 

Surface 

Inner  Surface 

Roughness 

- 

Medium 

Smooth 

Smooth 

Smooth 

Smooth 

Smooth 

Thermal  absorptance 

Fraction 

0.897 

0.900 

0.799 

0.780 

0.72  V  0.742b 

Solar  absorptance 

Fraction 

0.571 

0.9 

0.13 

0.13 

0.2 

Visible  absorptance 

0.571 

0.9 

0.13 

0.13 

0.2 

Thermal  resistance 

m2.K/W 

0.088 

0.46^= 

0.5V  0.0176b 

a  Surface  property  of  polystyrene  for  cooling  season  tests 
b  Surface  property  of  standard  floor  material  for  heating  season  tests 
c  Bondcote  specification  sheet  for  Extreme  250  insulating  fabric  with  Thinsulate 


3.1.3  Shade  fly  thermal  performance 

The  shade  fly  for  both  shelters  is  a  tan  colored  mesh  that  is  supported  6  - 
12  in.  above  the  shelter.  The  detached  shade  model  in  EnergyPlus  is  a  sim¬ 
ple  model  created  for  solid  or  partial  shading  of  direct  and  diffuse  solar 
radiation.  However,  it  is  thought  to  be  inadequate  for  the  shade  fly  on  the 
shelters  for  several  reasons: 

•  The  shading  model  does  not  include  transmittance  as  a  function  of  so¬ 
lar  angle  of  incidence. 

•  The  shading  model  does  not  include  the  partial  blocking  of  infrared 
heat  exchange  that  occurs  between  the  outer  fabric  of  the  shelter  and 
the  sky. 

•  The  shading  model  does  not  include  the  effect  that  solar  radiation 
striking  the  shading  fabric  will  have  on  the  infrared  heat  exchange  be¬ 
tween  the  outer  surface  of  the  shelter  and  the  shading  fabric.  The  shad¬ 
ing  fabric  will  be  warm  when  irradiated  by  the  sun,  which  will  reduce 
the  infrared  heat  loss  from  the  shelter  fabric  to  the  shading  fabric. 

•  The  shading  model  does  not  correctly  model  the  convective  thermal 
conditions  between  the  outer  surface  of  the  shelter  and  the  shading 
fabric. 

New  thermal  models  were  developed  for  the  shade  fly  on  each  shelter  in  an 
attempt  to  model  these  heat-transfer  paths.  These  models  are  described  in 
sections  3. 1.3.1  and  3.1.3.2,  and  are  implemented  in  the  EnergyPlus  Ener¬ 
gy  Management  System  (EMS)  code  that  sits  outside  of  the  EnergyPlus 
simulations. 
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3.1.3.1  Heat  exchange  through  shading  fabrics 

Figure  30  shows  the  thermal  network  used  to  calculate  convective  and 
radiative  heat  loss  from  the  outside  surface  of  a  shelter  when  a  fabric  is 
placed  as  a  shade  fly  over  the  shelter.  The  total  convective  and  radiative 
heat  loss  from  the  shelter  surface  in  the  “correct”  circuit  is  equal  to  the  to¬ 
tal  convective  and  radiative  heat  loss  in  the  “equivalent”  circuit.  At  each 
time  step  in  the  simulation,  the  “correct  circuit”  is  iteratively  solved  in  the 
EMS  code  and  then  used  to  calculate  values  for  the  convection  heat- 
transfer  coefficient,  emissivity,  and  amount  of  solar  radiation  reaching  the 
shelter’s  outer  surface  in  the  “equivalent  circuit.”  These  calculated  heat- 
transfer  coefficients  are  used  in  the  EnergyPlus  simulation. 


Figure  30.  Thermal  models  for  heat  transfer  through  the  shading  fabric. 

Shelter  Outer  Surface 


The  following  is  a  key  to  the  notation  in  Figure  30 

di  =  distance  between  edges  of  square  holes  in  shading  fabric 
da  =  length  of  side  of  one  square  hole  in  shading  fabric 
Tsi  =  temperature,  shelter  surface 

Tai  =  temperature,  air  gap  between  shelter  surface  and  shading 
fabric 
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Ta2  =  temperature,  outdoor  air 
Tsky  =  temperature,  effective  sky 
Tgnd  =  temperature,  ground  surface 
hci  =  convective  heat-transfer  coefficient,  shelter  surface  to  air  in 

gap 

hc2  =  convective  heat-transfer  coefficient,  air  in  gap  to  inside  surface 
of  shading  fabric 

hcinf  =  heat-transfer  coefficient,  gap  to  outdoor  air  (i.e.  infiltration 
into  gap) 

hc3  =  convective  heat-transfer  coefficient,  outside  surface  of  shading 
fabric  to  outdoor  air 

hceq  =  equivalent  convective  heat-transfer  coefficient,  outside  surface 
of  shading  fabric  to  outdoor  air 

hm  =  radiative  heat-transfer  coefficient,  shelter  surface  to  inside 
surface  of  shade 

hR2  =  radiative  heat-transfer  coefficient,  outside  shade  surface  to  sky 
hR3  =  radiative  heat-transfer  coefficient,  outside  shade  surface  to 
ground 

hR4  =  radiative  heat-transfer  coefficient,  shelter  surface  to  sky 
Hrs  =  radiative  heat-transfer  coefficient,  shelter  surface  to  ground 
Hrs  =  radiative  heat-transfer  coefficient,  outside  shade  surface  to 
outdoor  air 

hRy  =  radiative  heat-transfer  coefficient,  shelter  surface  to  outdoor 
air 

hReqi  =  equivalent  radiative  heat-transfer  coefficient,  shelter  surface  to 
outdoor  air 

hReq2  =  equivalent  radiative  heat-transfer  coefficient,  shelter  surface  to 
sky 

h-Reqs  =  equivalent  radiative  heat-transfer  coefficient,  shelter  surface  to 
ground 

fh  =  hole  fraction  of  shade  fabric 
a  =  solar  absorptance  of  shelter  surface 
Ic  =  global  solar  radiation  incident  on  shelter  surface 
T  =  true  overall  solar  transmittance  of  shading  fabric 
Teq  =  equivalent  solar  transmittance  of  shading  fabric 

The  radiative  heat-transfer  coefficients  are  defined  as: 


/i«  =  ae^FiTl  -  T^XT^  -K  T^) 


(9) 
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where 

(7  =  Stefan-Boltzmann  constant 
Si  =  emissivity  of  surface  i 
F  =  view  factor  of  surface  i  to  surface  2 
Ti  =  temperature  of  surface  1 
T2  =  temperature  of  surface  2 

The  total  rate  of  energy  leaving  the  node  Tsi  is  calculated  as  two  different 
values:  the  first  (Qi)  as  if  the  shading  fabric  has  no  holes,  the  second  (Q2) 
as  if  the  fabric  is  not  there  (i.e.,  “all  holes”).  The  total  rate  of  convective 
and  radiative  loss  from  Tsi  is  then  calculated  in  Equation  10  as  a  weighted 
average  of  the  two  values,  weighted  by  the  hole  fraction,  fh: 

Qtot  =  QiiX-fh)  +  Q2{fh)  (10) 

It  should  also  be  noted  that  in  our  analysis  we  have  assumed  a  high 
enough  value  of  hdnf  such  that  Tai  is  almost  equal  to  Ta2,  implying  a  very 
large  infiltration  rate  into  the  gap  between  the  shelter  surface  and  the 
shading  fabric. 

3. 1.3.2  Solar  transmittance  of  shade  fabric 

The  results  of  solar  transmittance  tests  on  the  shading  material  are  shown 
in  Figure  31  and  Figure  32  for  the  Utilis  and  AirBeam  shelters.  The  graphs 
show  the  strong  dependence  of  beam  solar  transmittance  on  the  incidence 
angle.  The  data  were  collected  by  simultaneously  measuring  the  beam  and 
diffuse  solar  radiation  on  a  clear  day  in  front  of  and  behind  a  sample  of  the 
shading  fabric  at  various  angles  relative  to  the  sun’s  position. 
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Figure  31.  Solar  transmittance  as  a  function  of 
beam  incidence  angle  for  the  Utilis  shading  fabric. 
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Figure  32.  Solar  transmittance  as  a  function  of 
beam  incidence  angle  for  the  AirBeam  shading  fabric. 
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A  new  model  was  created  to  approximate  the  effects  of  the  shade  fly  mesh 
that  takes  into  account  the  angular  dependence  of  the  solar  transmittance. 
This  algorithm  was  implemented  in  the  EnergyPlus  EMS  code  as  Equa¬ 
tions  11  and  12: 


“  fb^dn  fbd^dn  fd^d 

(11) 

^sh  ~  ^g,sh/^g 

(12) 

where 

Ig  sh  =  global  solar  radiation  on  the  shelter  surface  [W/m^] 

Idn  =  direct  normal  solar  radiation  on  the  shading  surface  [W/m^] 

Id  =  total  diffuse  (sky  and  ground)  solar  radiation  on  the  shading 
surface  [W/m^] 

Ig  =  global  solar  radiation  on  the  shading  surface  [W/m^] 

6i  =  incidence  angle  of  direct  normal  solar  radiation 
fd  =  fraction  of  transmitted  through  the  shade  as  beam  radiation 
fdd  =  fraction  of  Idn  transmitted  through  the  shade  as  diffuse 
radiation 

fd  =  fraction  of  Id  transmitted  through  the  shade 
Tsh  =  effective  transmittance  of  the  shade  fly 

Equation  ii  estimates  the  global  solar  radiation  that  is  transmitted 
through  the  shade  fly  mesh  and  is  incident  on  the  outer  surface  of  the  shel¬ 
ter.  Equation  12  is  the  effective  solar  transmittance  that  is  calculated  as  the 
ratio  of  the  global  solar  radiation  transmitted  through  the  shade  fly  to  the 
global  solar  radiation  on  the  shade  fly.  This  transmittance  is  calculated  in 
the  EnergyPlus  EMS  code  and  inserted  into  the  EnergyPlus  simulation  at 
each  time  step  during  daytime  hours. 

The  fractions  of  the  transmitted  solar  radiation  were  derived  from  the 
graphs  in  Eigure  31  and  Figure  32.  The  equations  for  these  correlations  are 
shown  in  Equations  13  -  15,  and  the  coefficients  are  shown  in  Table  5. 


fd  —  <^0 


(13) 


fd  =  &i(90  -  9i)  +  &2(90  -  OtY  -K  bs(.90  -  (14) 


fbd  =  cd90  -  dO  +  C2i90  -  di^ 


(15) 
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Table  5.  Coefficients  for  solar  fraction  correlations. 


Coefficient 

Utilis  (tan) 

HOT  (tan) 

ao 

0.45 

0.36 

bi 

3.104E-03 

5.747E-03 

b2 

-4.681E-06 

-6.257E-05 

bs 

-9.193E-08 

2.189E-07 

Cl 

6.759E-03 

4.824E-03 

C2 

-4.217E-05 

-2.954E-05 

3.1.4  Details  specific  to  model  calibration  study 

Test  data  were  stored  as  15  minute  averages  or  totals  for  all  measurement 
points.  All  simulations  were  run  on  a  15  minute  time  step.  Schedules  on  15 
minute  time  intervals  of  infiltration  and  indoor  air  temperature  were  cre¬ 
ated  from  the  measured  data  and  used  as  inputs  to  the  models.  The  time 
series  of  indoor  air  temperature  was  used  in  the  models  as  the  thermostat 
heating  set  point. 

The  model  of  the  floor  of  each  shelter  for  calibration  runs  was  designed  to 
ensure  that  the  modeled  total  heat  flow  through  the  floor  would  be  equal 
to  the  measured  heat  flow  (see  Equation  5).  The  approach  used  was  to 
model  a  fictitious,  massless  insulating  layer  of  arbitrary  thermal  resistance 
(0.5  m2-K/W)  on  the  floor  with  a  calculated  temperature  on  the  bottom 
side  of  the  insulating  layer  for  each  15  minute  time  step.  This  is  shown  in 
Equation  16: 


Tbot  —  T'f~  Qf^ 


(16) 


where 

Tbot  =  temperature  of  bottom  surface  of  fictitious  insulating  layer 
[°C] 

Tf  =  measured  floor  surface  temperature  [°C] 

Of  =  measured  floor  heat  flux  [W/m^] 

R  =  thermal  resistance  of  fictitious  layer  [m^-K/W] 

A  schedule  of  values  of  Tbot  for  each  15  minute  time  interval  during  the  test 
was  created  using  Equation  14.  Using  this  temperature  schedule  as  an  in¬ 
put,  the  model  then  calculated  the  resulting  top  surface  temperature  and 
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resulting  heat  flow  through  the  floor.  If  the  model  were  perfect,  the  mod¬ 
eled  top  surface  temperature  would  equal  the  measured  temperature  at 
each  simulation  time  step,  and  the  heat  flow  to  the  ground  would,  by  defi¬ 
nition,  match  the  measured  heat  flux.  Trying  to  match  the  modeled  floor 
surface  temperature  to  the  measured  temperature  was  part  of  the  model- 
calibration  process. 

During  the  heating-season  tests  the  indoor  air  temperature  was  main¬ 
tained  by  controlling  convective  electric  heaters.  To  emulate  this  type  of 
heating  in  the  model  for  the  purposes  of  model  calibration,  the  “Base¬ 
board  :RadiantConvective”  model  available  in  EnergyPlus  was  used,  with  a 
radiant  fraction  of  o.o. 

During  the  cooling-season  tests,  a  5  ton  ECU  was  used  to  maintain  the 
cooling  set  point.  To  emulate  this  type  of  cooling  device  in  the  model,  the 
“Packaged  Terminal  Air  Conditioner”  (PTAC)  model  was  used  in 
EnergyPlus  to  represent  the  ECU.  The  model  merely  had  to  maintain  the 
air  temperature  of  the  interior  of  the  structure  by  removing  sensible  ener¬ 
gy  from  the  air,  and  the  magnitude  of  this  energy  removal  was  compared 
to  the  measured  value  (Equation  8).  Therefore,  the  PTAC’s  modeled  effi¬ 
ciency  was  not  relevant  to  calibrating  the  thermal  performance  of  the  shel¬ 
ter  itself. 

Custom  fifteen-minute-interval  weather  files  were  generated  with  the 
EnergyPlus  Weather  Converter  program  using  one-minute  averaged 
weather  station  data  collected  during  each  field  test.  The  following  meas¬ 
urements  were  included: 

•  Outdoor  dry  bulb  temperature 

•  Relative  humidity 

•  Wind  speed 

•  Wind  direction 

•  Sky  radiation 

•  Global  horizontal  solar  radiation 

•  Direct  normal  solar  radiation 
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3.2  Calibration  study 

3.2.1  Approach 

The  general  philosophy  of  model  calibration  is  to  limit  the  calibration  to 
the  adjustment  of  parameters  that  cannot  be  known  with  enough  certainty 
to  be  incorporated  into  the  model.  For  example,  the  solar  absorptance  of 
the  outer  fabric  of  the  shelter  was  measured  in  a  laboratory  at  NREL.  By 
measuring  difficult-to-predict  heat  flows  (infiltration,  ground  heat  flux) 
the  calibration  process  was  made  simpler:  the  heat  flows  not  measured 
were  the  heat  flows  through  the  walls  and  roof  of  the  shelter,  so  that  only 
those  parameters  that  affect  the  heat  flow  through  the  roof  and  walls  need 
to  be  adjusted  to  match  model  to  measurement.  Among  those  parameters, 
many  are  constants  that  were  measured  in  the  laboratory  (i.e.,  fabric  prop¬ 
erties  listed  in  Table  3  and  Table  4). 

Figure  33  is  a  simplified  diagram  of  the  thermal  network  in  the  Utilis 
model.  Resistances  to  heat  flow  are  shown  as  resistor  symbols  as  used  in 
electronics  circuit  diagrams.  Resistances  between  a  surface  temperature 
and  an  air  temperature  represent  resistance  to  convective  heat  flow,  gov¬ 
erned  by  convective  heat-transfer  coefficients.  Resistances  between  two 
surface  temperatures  represent  resistance  to  infrared  heat  flow,  governed 
by  surface  emissivities.  The  emissivities  of  the  surfaces  were  measured  in 
the  laboratory,  leaving  the  most  significant  uncertain  heat-flow  parame¬ 
ters  in  the  model  to  be  the  convection  coefficients  on  each  fabric  surface. 


Figure  33.  Simplified  thermal  network  showing 
resistances  to  heat  flow  through  various  paths. 
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In  typical  building  simulations,  the  exact  magnitudes  of  convective  heat- 
transfer  coefficients  are  not  very  important  because  the  dominant  re¬ 
sistance  to  heat  flow  is  the  thermal  conductivity  of  the  walls.  In  the  case  of 
fabric  structures,  however,  the  resistance  to  thermal  conduction  through 
the  fabric  itself  is  low  so  the  convective  and  radiative  heat  transfers  are  the 
dominant  mechanisms  for  heat  flow.  However,  the  correlations  for  convec¬ 
tive  heat-transfer  coefficients  in  EnergyPlus  were  developed  for  rigid  sur¬ 
faces  found  in  buildings,  and  are  not  appropriate  for  moving  fabric 
surfaces.  A  study  by  Zaki,  Al-Turki,  and  Fatani  (1993)  estimated  that  the 
convective  heat-transfer  coefficients  on  taut  fabric  surfaces  could  be  ap¬ 
proximately  two  to  three  times  the  order  of  magnitude  of  those  for 
smooth,  rigid  surfaces.  The  surfaces  in  these  military  shelters  are  relatively 
loose  and  are  moved  freely  by  an  ECU  fan  and  wind,  potentially  causing 
even  higher  convection  heat-transfer  coefficients.  In  addition,  the  shapes 
of  the  shelters  are  not  as  precisely  aligned  in  the  EFOB-L  as  the  drawings 
here  may  imply.  The  air  gaps  between  the  inner  and  outer  fabrics  are 
shaped  by  draped  fabrics  enclosing  irregular  air  pockets.  The  natural  con¬ 
vection  patterns  that  develop  in  the  air  gaps  are  most  likely  not  well- 
behaved  as  the  convection  patterns  upon  which  the  standard  correlations 
are  based. 

Developing  new  correlations  for  loose  fabric  surfaces  was  beyond  the 
scope  of  this  project.  Therefore,  we  selected  the  best  correlations  from  the 
EnergyPlus  library  for  different  surface  types,  surface  orientations,  and 
direction  of  heat  flow  and  applied  scaling  factors  for  each  surface  type  that 
could  be  adjusted  to  calibrate  the  models  with  the  measured  energy  data. 

Comparing  model  to  measurement  requires  an  output  or  outputs  from  the 
simulation  that  can  be  compared  to  measured  values.  In  this  case  the  prin¬ 
cipal  measured  quantity  for  comparison  was  the  total  heating  or  cooling 
needed  to  maintain  the  shelter’s  indoor  air  temperature.  Secondary  meas¬ 
ured  quantities  included  the  surface  temperatures  of  the  fabrics.  These 
secondary  comparison  points  are  useful  in  deciding  whether  or  not  the 
correct  convection  heat-transfer  coefficients  are  being  used  on  each  sur¬ 
face.  If  only  the  total  heating  or  cooling  energy  is  used  as  a  basis  for  com¬ 
parison  between  model  and  measurement  there  would  likely  be  more  than 
one  combination  of  adjustments  to  the  many  convection  coefficients  that 
would  yield  similar  agreement  of  total  energy  between  model  and  meas¬ 
urement  for  the  test  period.  Some  of  these  combinations  would  not  neces¬ 
sarily  yield  good  agreement  between  model  and  measurement  under  a 
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different  set  of  test  conditions.  Since  the  goal  in  calibrating  the  model  is  to 
gain  confidence  in  a  model’s  ability  to  correctly  predict  performance  over  a 
wider  range  of  conditions  than  were  encountered  during  field  measure¬ 
ments,  it  was  important  to  have  the  secondary  comparisons  to  be  sure  all 
parts  of  the  model  are  behaving  correctly. 

The  energy  models  are  made  up  of  many  surfaces  as  shown  previously  in 
Figure  28  and  Figure  29,  and  each  has  an  associated  convection  coeffi¬ 
cient.  To  simplify  the  calibration  process,  we  broke  these  “faces”  into  five 
categories: 

•  Floor  surface 

•  Inside  face  of  inside  liner 

•  Outside  face  of  inside  liner  (inside  the  air  gap) 

•  Inside  face  of  outer  shell  (inside  air  gap) 

•  Outside  face  of  outer  shell 

Five  custom  scaling  factors,  each  applying  to  the  convection  coefficients 
associated  with  each  surface  category  were  developed  by  calibrating  the 
models  for  each  shelter  to  the  measured  data.  EnergyPlus  does  not  allow 
the  convection  heat-transfer  coefficients  to  be  adjusted  or  scaled;  there¬ 
fore,  the  convection  heat-transfer  coefficients  are  calculated  outside  of  the 
EnergyPlus  simulation  in  the  EnergyPlus  EMS  code.  Then  the  calculated 
heat-transfer  coefficients  with  the  appropriate  scaling  factors  are  inserted 
into  the  EnergyPlus  simulation  for  each  surface  at  each  time  step. 

The  scaling  factors  were  determined  during  the  model-calibration  process 
with  the  measured  field  data  in  a  custom  optimization  routine  built  in  Mi¬ 
crosoft  Excel  Visual  Basic  for  Applications  (VBA).  The  objective  function 
was  the  root  mean  square  (RMS)  errors  between  modeled  and  measured 
heating  load  and  the  modeled  and  measured  surface  temperatures.  The 
optimization  routines  were  run  to  minimize  the  RMS  errors  to  determine 
the  values  of  the  five  independent  variables  (i.e.,  the  scaling  factors).  These 
scaling  factors  were  then  used  for  the  annual  simulations. 

3.2.2  Results 

The  models  were  calibrated  on  the  measured  energy  flows  of  the  ECU  and 
the  ground  heat  flux  by  adjusting  the  convection  heat-transfer  coefficient 
scaling  factors  across  the  March  and  June  calibration  periods.  The 
AirBeam  model  was  calibrated  in  the  baseline  condition  because  that  is  the 
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configuration  that  was  field  tested.  The  Utilis  shelter  was  field  monitored 
with  a  shade  fly  and  a  Thinsulate  liner  with  an  added  radiant  barrier,  so 
the  calibration  energy  model  included  these  features.  We  attempted  to  cal¬ 
ibrate  the  Utilis  energy  model  with  the  newly  derived  models  described  in 
section  3.1.3.1  and  section  3. 1.3.2;  however,  the  calibration  results  were 
not  physically  sensible  and  resulted  in  high  errors  and  unrealistic  scaling 
factors.  We  were  unable  to  determine  the  source  of  these  errors  and  in¬ 
stead  used  the  EnergyPlus  detached  shading  model  with  a  solar  transmit¬ 
tance  of  0.44. 

Table  6  shows  the  average  ECU  energy  and  the  ground  heat  flux  with  the 
model  errors.  Table  7  shows  the  scaling  factors  that  are  applied  to  the  sur¬ 
face  convection  heat-transfer  coefficients  calculated  during  the  annual 
simulations. 


Table  6.  Model  calibration  resuits. 


Shelter 

Total  ECU  Energy 

Ground  Heat  Flux 

March 

June 

March 

June 

Avg(kW) 

RMS  Error 

Avg(kW) 

RMS  Error 

Avg(kW) 

RMS  Error 

Avg(kW) 

RMS  Error 

Air  Beam 

5.37 

9.2% 

6.36 

7.9% 

0.67 

46.9% 

-0.07 

34.9% 

Utilis 

2.89 

15.6% 

2.59 

10.6% 

0.93 

15.8% 

0.20 

10.2% 

Tabie  7.  Convection  heat-transfer  coefficient  scaiing  factors 


Shelter 

Inner  Shell 

Outer  Shell 

Floor 

Inner  Face 

Outer  Face 

Inner  Face 

Outer  Face 

Inner  Face 

Air  Beam 

4.58 

13.62 

7.01 

0.50 

1.25 

Utilis 

5.14 

6.37 

11.19 

0.50 

0.50 

Use  of  the  calibrated  scaling  factors  is  not  ideal,  but  it  does  allow  the  ener¬ 
gy  simulation  to  capture  some  of  the  dynamic  effects  of  the  convection  co¬ 
efficients  with  changing  conditions  and  provides  an  adjustment  for  the 
fabric  surface.  These  scaling  factors  are  derived  for  these  specific  shelter 
configurations  and  conditions,  so  they  may  not  translate  very  well  to  other 
applications. 


ERDC/CERLTR-15-13 


45 


4  Simulations  and  Analyses 

The  calibrated  models  were  used  to  assess  the  performance  of  alternative 
envelope  configurations  in  four  locations  to  determine  potential  energy 
savings.  These  locations  included  the  EFOB-L  at  Champaign,  IL;  Fort 
Devens,  MA;  Ali  Al  Salem,  Kuwait;  and  Panama  City,  FL. 

4.1  Annual  model  descriptions 

Adjustments  were  made  to  the  calibrated  energy  models  for  the  annual 
simulations.  Changes  to  the  calibrated  models  are  explained  below. 

4.1.1  Ground  heat  transfer 

The  floor  was  modeled  as  i  meter  of  soil  with  a  half  cm  air  gap  and  the 
floor  liner.  However,  the  combination  of  the  thin  floor  liner  and  the  ther¬ 
mally  massive  soil  caused  the  simulation  to  be  unstable  and  stop  prema¬ 
turely,  and  the  liner  was  eliminated  from  the  construction.  The  annual 
model  uses  average  monthly  ground  temperatures  at  i  m  depth,  which 
were  interpolated  from  the  monthly  ground  temperatures  at  0.5  m,  2  m, 
and  4  m  from  the  FnergyPlus  statistical  weather  information  file  for  each 
location.  The  floor  material  properties  are  listed  in  Table  8. 


Table  8.  Floor  material  properties. 


Property 

Units 

1  meter  soil 

Half-inch  air  gap 

Roughness 

- 

Smooth 

Smooth 

Thickness 

m 

1 

0.005 

Conductivity 

W/m-K 

0.88 

Density 

Kg/m3 

1200 

Specific  heat 

J/kg-K 

2192 

Thermal  absorptance 

Fraction 

0.9 

0.9 

Solar  absorptance 

Fraction 

0.5 

0.45 

Visible  absorptance 

Fraction 

0.5 

0.6 

Thermal  resistance 

m2.K/W 

- 

0.47 
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4.1.2  Infiltration 

Infiltration  was  modeled  in  Energy  Plus  using  the  Zoneinfiltration: 
EffectiveLeakageArea  object  following  Equation  i  and  the  equivalent  leak¬ 
age  areas  (ELA)  measured  during  the  summer  monitoring  period  using  a 
blower-door  analysis  as  shown  in  Table  9.  The  stack  and  wind  coefficients 
were  taken  from  the  EnergyPlus  Input/Output  Reference  for  a  one-story 
house  with  no  obstructions  or  shielding  of  the  wind  (DOE  2014b).  We  do 
not  know  the  application  of  the  shelters  in  the  field,  but  we  selected  the 
no-shielding  case  because  the  shelters  are  sensitive  to  wind  speed  and  it  is 
the  most  conservative  scenario.  The  stack  and  wind  coefficients  for  Equa¬ 
tion  1  are  0.000145  and  0.000319,  respectively. 


Table  9.  Infiltration  ELAand  leakage  rates. 


Shelter 

ELA  (cm2) 

ELA  (in2) 

cfm/ft2  @4Pa 

cfm/ft2  @75Pa 

Air  Beam 

466 

41 

0.13 

0.70 

Util  is 

264 

72 

0.18 

0.98 

4.1.3  Thermostat  set  points  and  control 

The  assigned  heating  and  cooling  set  points  were  20  °C  (68  °E)  and  26.7 
°C  (80  °F),  respectively.  The  ECU  model  allowed  for  heating  or  cooling 
whenever  necessary. 

4.1.4  Internal  loads 

It  was  assumed  that  the  shelters  were  command  and  control  centers  with 
continuous  operation.  It  was  assumed  that  the  internal  electrical  equip¬ 
ment  loads  were  5  kW  and  that  there  were  20  occupants. 

4.2  Alternative  envelope  configurations 

The  calibrated  models  were  used  to  assess  the  performance  of  alternative 
envelope  configurations  to  determine  the  impact  on  ECU  energy  consump¬ 
tion.  The  alternative  envelope  options  are  listed  in  Table  10  and  the  mate¬ 
rial  properties  are  listed  in  Table  11. 
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Table  10.  Alternative  envelope  options. 


Envelope  Options 

Materials 

Shelter 

Baseline 

Inner  liner 

AirBeam;  Utilis 

Option  1 

Baseline  +  radiant  barrier 

AirBeam;  Utilis 

Option  2 

Baseline  +  shade  fly 

AirBeam;  Utilis 

Option  3 

Thinsulate  liner  * 

AirBeam;  Utilis 

Option  4 

Thinsulate  liner  and  radiant  barrier 

AirBeam;  Utilis 

Option  5 

Baseline  +  PCM 

AirBeam 

*  The  Thinsulate  liner  replaces  the  baseline  inner  liner. 


Tabie  11.  Enveiope  material  properties. 


Property 

Units 

Radiant 

Barrier 

Thinsulate 

PCM 

Shade 

Roughness 

- 

Medium 

Smooth 

Smooth 

Smooth 

- 

Thickness 

m 

- 

- 

0.01 

- 

Conductivity 

W/m-K 

- 

- 

0.2 

- 

Density 

K^m3 

- 

- 

850 

- 

Specific  heat 

J/kg-K 

- 

- 

2100 

- 

Thermal  absorptance 

Fraction 

0.05 

0.8 

0.9 

- 

Solar  absorptance 

Fraction 

0.05 

0.3 

0.7 

- 

Visible  absorptance 

Fraction 

0.05 

0.3 

0.7 

- 

Thermal  resistance 

m2.K/W 

0.03 

0.46 

- 

- 

Solar  Transmittance 

Fraction 

- 

- 

- 

0.42 

The  baseline  construction  for  both  shelters  consists  of  an  inner  liner  made 
of  a  white  polymer  canvas  and  an  outer  shell  made  of  a  slightly  heavier 
polymer  canvas.  The  color  of  the  outer  shell  is  tan  on  the  outside  and  black 
on  the  inside.  The  radiant  barrier  is  a  very  thin,  highly  reflective  layer  that 
lays  on  top  of  the  inner  liner.  The  Thinsulate  liner  is  sandwiched  between 
two  layers  of  the  inner  liner  material  with  a  very  thin  layer  of  a  radiant 
barrier  material. 

A  simulation  using  the  AirBeam  shelter  model  was  also  run  with  a  phase- 
change  material  (PCM)  liner  called  BioPCMat™  (Phase  Change  Energy 
Solutions,  Asheboro,  NC),  a  sheet-type  product  that  incorporates  pockets 
of  a  nontoxic,  biodegradable  soy-based  PCM  (BioPCM  2013).  Stock  shel¬ 
ters  of  the  type  modeled  have  near- zero  thermal  mass  except  for  contact 
with  the  ground,  so  the  addition  of  the  PCM  liner  was  expected  to  have  a 
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larger  impact  on  the  shelter  performance  than  it  would  in  a  conventional 
building.  The  thermodynamic  and  enthalpy  properties  (for  cooling  mode) 
are  listed  in  Table  ii  (above)  and  Table  12,  respectively.  The  temperature- 
enthalpy  curve  was  derived  from  experimental  measurements  of  tempera¬ 
ture  and  specific  heat  taken  for  a  concurrent  PCM  study  by  ERDC-CERL 
researchers  (results  being  prepared  for  publication  in  separate  report.  A 
thin  blanket  of  this  material  was  applied  to  the  outside  of  the  inner  liner. 
The  PCM  blanket  has  pockets  of  PCM  material,  and  it  was  assumed  that 
the  weight  of  PCM  material  was  0.53  Ib/ft^  (2.59  kg/m^).  The  shelter  with 
PCM  liner  was  simulated  in  the  eight  U.S.  climate  zones  listed  in 

Table  13. 


Table  12.  Temperature  vs  enthalpy  of  commercial  PCM. 


Temperature  (®C) 

Enthalpy  (J/kg) 

10.01 

0.0142 

15.18 

10.398 

21.09 

107.929 

21.77 

142.430 

22.05 

152.413 

23.04 

154.882 

26.61 

160.025 

Table  13.  U.S.  climate  zones  for  PCM  modeling  In  AIrBeam  shelter. 


Climate  Zone 

Climate  Description 

Representative  City 

lA 

Hot,  humid 

Miami,  FL 

2B 

Hot,  dry 

Phoenix,  AZ 

3B 

Warm,  dry 

Los  Angeles,  CA 

4B 

Mild,  dry 

Albuquerque,  NM 

5B 

Cold,  dry 

Boulder,  CO 

6B 

Cold,  dry 

Helena,  MT 

7 

Very  cold 

Duluth,  MN 

8 

Extremely  cold 

Fairbanks,  AK 

Locations 

Simulations  of  both  shelters  were  run  using  weather  data  for  Champaign, 
IL;  Fort  Devens,  MA;  and  Ali  A1  Salem,  Kuwait.  EnergyPlus  weather  files 
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derived  from  National  Solar  Radiation  Data  Base  TMY3  (Typical  Meteoro¬ 
logical  Year  3)  data  were  used  for  each  location.  In  addition,  as  noted 
above,  the  AirBeam  shelter  with  the  PCM  liner  was  simulated  in  eight  U.S. 
locations  representing  each  major  climate  zone  (shown  in  Table  13). 

4.4  Performance  results 

In  all  cases,  the  annual  energy  consumption  for  lighting  and  electrical 
equipment  was  43,800  kWh,  which  is  roughly  two  thirds  of  the  total  ener¬ 
gy  consumption.  The  lighting  and  equipment  along  with  the  occupants 
represents  a  significant  internal  load  that  has  to  be  removed  from  the  shel¬ 
ter  to  maintain  the  cooling  set  point.  The  load  is  so  high  that  it  results  in 
very  high  cooling  loads  even  in  the  cold  months  when  the  shelter  is  closed 
up  with  normal  amount  of  infiltration.  A  much  more  efficient  way  to  oper¬ 
ate  the  shelters  in  this  case  is  to  open  the  doors  and  windows  when  it  is 
cool  outside.  Natural  ventilation  was  included  in  the  EnergyPlus  models  at 
5  ACH  to  simulate  opening  the  doors  and  windows  when  the  outside  air 
temperature  was  between  o  °C  and  16.5  °C  (32  °F  and  62  °F).  The  natural 
ventilation  is  not  allowed  to  cool  the  shelter  below  the  heating  set  point  of 
20  °C.  All  of  the  simulations,  including  the  baseline,  applied  natural  venti¬ 
lation  to  provide  free  cooling  when  it  was  available.  Without  the  natural 
ventilation,  the  liners  result  in  higher  ECU  energy  consumption  than  the 
baseline  because  of  the  increased  cooling  loads  from  the  trapped  heat. 

4.4.1  AirBeam 

The  annual  performance  for  the  AirBeam  shelter  with  the  alternative  enve¬ 
lope  configurations  are  presented  in  Table  14,  Table  15,  and  Table  16,  and 
shown  graphically  in  Figure  34.  The  monthly  ECU  energy  totals  for  the 
baseline  and  alternative  (radiant  barrier)  scenarios  at  Champaign,  IL,  are 
shown  in  Figure  35  through  Figure  38. 


Table  14.  AirBeam  results  for  Champaign,  IL. 


Energy  Consumption 
(kWh) 

Baseline 

Radiant 

Barrier 

Th  insulate 
Liner 

PCM  Liner 

Shade  Fly 

Heating 

11,017 

2,914 

178 

2,253 

10,800 

Cooling 

11,681 

12,442 

13,681 

12,589 

10,611 

Fans 

5,022 

5,022 

5,022 

5,022 

5,022 

Total  ECU 

27,719 

20,377.8 

18,881 

19,864 

26,436 

ECU  Energy  Saved 

7,341.7 

8,839 

7,856 

1,283 

ECU  Percent  Savings 

26.5% 

31.9% 

28.3% 

4.6% 
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Table  15.  AirBeam  results  for  Fort  Devens,  MA. 


Energy  Consumption 
(kWh) 

Baseline 

Radiant 

Barrier 

Thinsulate 

Liner 

PCM  Liner 

Shade  Fly 

Heating 

13,511 

3,211 

142 

2,425 

13,381 

Cooling 

7,094 

8,306 

9,778 

8,519 

6,133 

Fans 

5,022 

5,022 

5,022 

5,022 

5,022 

Total  ECU 

25,631 

16,538.9 

14,942 

15,964 

24,536 

ECU  Energy  Saved 

9,091.7 

10,689 

9,667 

1,094 

ECU  Percent  Savings 

35.5% 

41.7% 

37.7% 

4.3% 

Table  16.  AirBeam  results  for  AN  Al  Salem,  Kuwait. 


Energy  Consumption 
(kWh) 

Baseline 

Radiant 

Barrier 

Thinsulate 

Liner 

PCM  Liner 

Shade  Fly 

Heating 

56 

- 

- 

- 

31 

Cooling 

44,831 

40,625 

40,047 

40,275 

43,244 

Fans 

5,022 

5,022 

5,022 

5,022 

5,022 

Total  ECU 

49,908 

45,650.0 

45,069 

45,297 

48,297 

ECU  Energy  Saved 

4,258.3 

4,839 

4,611 

1,611 

ECU  Percent  Savings 

8.5% 

9.7% 

9.2% 

3.2% 

Figure  34.  Annuai  ECU  energy  consumption  and  savings  for  the  AirBeam  sheiter. 
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Figure  35.  Monthly  ECU  energy  performance  for  the  AirBeam  shelter  in  Champaign,  IL 
for  the  baseline  (left  bar)  and  radiant  barrier  (right  bar)  cases. 


Figure  36.  Monthly  ECU  energy  performance  for  the  AirBeam  shelter  in  Champaign,  IL 
for  the  baseline  (left  bar)  and  Thinsulate  (right  bar)  liner  cases. 
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Figure  37.  Monthly  ECU  energy  performance  for  the  AirBeam  shelter  in  Champaign,  IL 
for  the  baseline  (left  bar)  and  radiant  barrier  and  Thinsulate  liner  (right  bar)  cases. 
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Figure  38.  Monthly  ECU  energy  performance  for  the  AirBeam  shelter  for  in 
Champaign,  IL  the  baseline  (left  bar)  and  the  shade  fly  (right  bar)  cases. 


The  Thinsulate  liner  provides  the  best  performance  in  every  climate.  How¬ 
ever,  the  radiant  barrier  and  PCM  liner  provide  similar  savings.  The  mesh 
shade  fly  provided  the  smallest  savings  in  all  climates.  We  were  not  able  to 
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test  the  baseline  shelter  configuration  with  and  without  the  shade  fly,  and 
therefore  were  not  able  to  validate  the  model  of  the  shade  fly. 

The  mesh  shade  fly  has  multiple  impacts  on  the  shelter  in  that  it  blocks 
some  of  the  incoming  solar  radiation,  but  it  also  creates  a  layer  of  warm  air 
near  the  outer  surface  and  blocks  nighttime  radiation  exchange  with  the 
sky.  During  the  cooling  season,  the  shade  fly  blocks  solar  gain  but  it  can 
also  cause  higher  convective  heat  gains  by  trapping  warm  air  near  the  sur¬ 
face  and  reduce  heat  loss  on  clear  nights  by  blocking  the  long- wave  radia¬ 
tion  exchange.  In  the  heating  season,  the  shade  fly  reduces  heat  load  by 
acting  as  another  barrier  to  convective  heat  loss  to  the  cold  ambient  air 
and  partially  blocks  the  long-wave  radiation  exchange  with  clear  night 
skies,  but  it  also  blocks  beneficial  solar  gains  on  sunny  days.  The  optimal 
performance  for  the  shade  fly  would  be  to  only  install  it  during  cooling 
season  and  remove  it  during  the  heating  season,  but  that  scenario  was  not 
simulated  in  this  project. 

The  monthly  results  for  Champaign,  IL  show  the  variation  in  the  impact 
from  the  different  liners  on  the  heating  and  cooling  loads.  The  Thinsulate 
liner  nearly  eliminates  the  need  for  heating  because  most  heating  is  pro¬ 
vided  by  the  internal  loads.  The  cooling  loads  increase  in  every  case  except 
for  the  shade  fly,  which  reduces  the  cooling  load  compared  to  the  baseline. 

It  should  be  noted  at  this  point  that  the  results  for  a  shelter  used  for  sleep¬ 
ing  quarters  would  differ  greatly  because  the  internal  gains  would  be  dras¬ 
tically  reduced.  The  cooling  loads  would  be  much  lower  and  the  heating 
loads  would  be  much  higher.  The  savings  of  each  of  the  liners  and  shade 
fly  would  be  different  as  well.  The  Thinsulate  liner  with  radiant  barrier 
would  probably  be  the  best  performer,  but  that  is  not  certain  without  addi¬ 
tional  energy  modeling  or  testing. 

The  performance  results  for  the  AirBeam  with  the  PCM  liner  are  shown  in 
Table  17  and  Figure  39  for  the  eight  supplemental  locations.  The  impact  of 
the  PCM  liner  on  annual  ECU  energy  consumption  is  similar  to  the 
Thinsulate  liner  (see  Figure  34).  The  PCM  liner  slightly  increases  the  cool¬ 
ing  load  in  most  climates  but  has  a  significant  impact  on  heating  energy. 

In  this  application,  the  PCM  liner  acts  as  an  insulating  layer;  it  absorbs  the 
heat  from  the  internal  gains  and  slowly  releases  the  energy  to  the  air  gap 
between  the  inner  liner  and  the  outer  shell.  The  results  for  the  PCM  liner 
should  be  taken  to  be  preliminary  because  the  PCM  model  in  EnergyPlus 
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has  not  been  rigorously  validated  at  this  point,  and  the  PCM  liner  was  not 
part  of  the  field  monitoring  in  this  project. 


Table  17.  AirBeam  ECU  energy  consumption. 


lA  Miami 

2B  Phoenix 

3B  Los  Angeles 

4B  Albuquerque 

5B  Denver 

6B  Helena 

7  Duluth 

8  Fairbanks 

Baseline 

40,131 

37,942 

17,303 

23,906 

22,139 

22,011 

30,647 

35,722 

PCM  Liner 

39,522 

35,375 

21,006 

20,322 

16,486 

14,150 

16,903 

19,522 

Savings 

608 

2,567 

-3,703 

3,583 

5,653 

7,861 

13,744 

16,200 

%  Savings 

1.5% 

6.8% 

-21.4% 

15.0% 

25.5% 

35.7% 

44.8% 

45.3% 

Figure  39.  Annuai  ECU  energy  consumption  and  savings 
for  the  AirBeam  sheiter  with  the  PCM  iiner. 


4.4.2  Utilis 

The  results  of  for  the  Utilis  shelter  with  the  alternative  envelope  configura¬ 
tions  are  presented  in  Table  i8  through  Table  20,  and  shown  graphically  in 
Figure  40.  The  radiant  liner  and  the  Thinsulate  liner  provide  nearly  the 
same  annual  energy  performance  for  the  two  heating-dominated  climates, 
but  for  slightly  different  reasons.  The  Thinsulate  liner  provides  better  re¬ 
duction  of  the  heating  loads,  but  the  radiant  liner  provides  the  best  reduc¬ 
tion  in  cooling  loads.  Combining  the  radiant  barrier  to  the  Thinsulate  liner 
eliminates  the  heating  load,  but  increases  the  cooling  load  and  provides 
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smaller  savings.  This  occurs  because  it  traps  more  of  the  internal  gains 
thus  requiring  more  cooling.  The  Thinsulate  liner  with  radiant  barrier  is 
the  highest  performer  in  the  very  hot  climate  because  it  is  the  most  effec¬ 
tive  at  blocking  the  solar  and  conductive  heat  gains.  The  mesh  shade  fly 
provided  the  smallest  savings  in  all  climates  except  for  Panama  City,  FL. 
The  shade  fly  works  well  for  the  moderate  climates  but  not  for  the  extreme 
climates.  The  performance  would  probably  be  improved  if  it  were  only  in¬ 
stalled  in  the  cooling  months. 


Table  18.  Utilis  results  for  Champaign,  IL 


Energy  Consumption 
(kWh) 

Baseline 

Radiant 

Barrier 

Thinsulate 

Liner 

Radiant  & 
Thinsulate 

Shade  Fly 

Heating 

8,003 

700 

56 

- 

7,539 

Cooling 

13,653 

14,444 

15,117 

17,964 

11,825 

Fans 

7,825 

7,825 

7,825 

7,825 

7,825 

Total  ECU 

29,481 

22,969.4 

22,997 

25,792 

27,192 

ECU  Energy  Saved 

6,511.1 

6,483 

3,689 

2,289 

ECU  Percent  Savings 

22.1% 

22.0% 

12.5% 

7.8% 

Table  19.  Utilis  resuits  for  Fort  Devens,  MA. 


Energy  Consumption 
(kWh) 

Baseline 

Radiant 

Barrier 

Thinsulate 

Liner 

Radiant  & 
Thinsulate 

Shade  Fly 

Heating 

9,425 

672 

67 

- 

8,919 

Cooling 

9,117 

10,328 

11,289 

14,781 

7,331 

Fans 

7,825 

7,825 

7,825 

7,825 

7,825 

Total  ECU 

26,367 

18,827.8 

19,183 

22,608 

24,078 

ECU  Energy  Saved 

7,538.9 

7,183 

3,758 

2,289 

ECU  Percent  Savings 

28.6% 

27.2% 

14.3% 

8.7% 

Tabie  20.  Utiiis  resuits  for  Aii  Ai  Saiem,  Kuwait. 


Energy  Consumption 
(kWh) 

Baseline 

Radiant 

Barrier 

Thinsulate 

Liner 

Radiant  & 
Thinsulate 

Shade  Fly 

Heating 

42 

- 

- 

- 

17 

Cooling 

48,128 

43,361 

42,128 

40,836 

45,264 

Fans 

7,825 

7,825 

7,825 

7,825 

7,825 

Total  ECU 

55,997 

51,186.1 

49,953 

48,664 

53,106 

ECU  Energy  Saved 

4,811.1 

6,044 

7,333 

2,892 

ECU  Percent  Savings 

8.6% 

10.8% 

13.1% 

5.2% 
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Figure  40.  Annual  ECU  energy  consumption  and  savings  for  the  Utilis  shelter. 
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5  Conclusions 

5.1  Conclusions 

This  project  demonstrated  that  whole-building  energy  monitoring  and 
modeling  techniques  can,  with  minor  modifications,  be  used  to  develop 
calibrated  simulation  models  of  energy  performance  for  soft-walled  mili¬ 
tary  shelters.  Some  uncertainties  in  the  demonstrated  models  impose  limi¬ 
tations  on  their  applicability,  but  the  usefulness  of  these  calibrated  energy 
models  has  been  demonstrated  for  purposes  of  showing  the  relative  energy 
performance  of  different  shelter  configurations,  operations,  and  locations. 

The  results  of  the  energy  modeling  show  that  for  a  command  operations 
type  of  shelter  with  very  high  internal  loads,  the  radiant  barrier  or  the 
Thinsulate  liner  alone  provide  the  best  performance  in  cold  climates;  and 
the  Thinsulate  liner  with  the  radiant  barrier  provides  the  best  performance 
in  the  hot  climates.  The  shade  fly  provided  the  least  energy  savings  in  all 
climate  simulations  except  for  Panama  City,  FL,  where  it  provided  the 
highest  savings.  The  results  would  probably  be  improved  if  the  shade  fly 
were  installed  only  during  hot  months  and  removed  for  the  cooler  months. 

The  simulations  predict  that  the  phase-change  material  (PCM)  liner  will 
provide  significant  savings  in  the  cold  climates,  although  not  as  much  as 
the  Thinsulate  liner.  The  results  for  the  PCM  liner  should  be  considered 
preliminary  because  we  did  not  take  field  measurements  of  the  liner’s  per¬ 
formance,  and  at  this  point  there  has  been  very  little  validation  of  the  PCM 
model  in  EnergyPlus. 

The  results  presented  in  this  report  are  highly  dependent  on  the  assump¬ 
tions  for  shelter  operation,  including  the  internal  loads,  and  how  tightly 
the  shelter  is  sealed  against  infiltration.  The  results  depend  greatly  on  the 
use  of  natural  ventilation  instead  of  ECU  cooling  when  it  is  advantageous. 
Actual  results  in  the  field  may  differ  significantly  based  on  the  actual  shel¬ 
ter  operation,  ECU  performance,  and  local  weather  conditions.  However, 
we  are  confident  in  the  relative  performance  and  the  trends  shown  in  the 
simulation  results  for  the  different  shelter  configurations. 
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5.2  Recommendations 

It  is  recommended  that  a  tactical  soft  wall  shelter  energy  “virtual  test  bed” 
be  developed  for  DoD  soft  wall  shelters  with  Energy  Plus  modeling,  simu¬ 
lation  and  validation  to  reduce  the  number  of  field  evaluation  tests,  since  a 
validated  model  can  predict  the  shelter  performance  with  various  linings 
(with  known  properties)  in  all  climate  zones  to  mitigate  energy  losses. 

New  shelter  linings  and  shelter  designs  based  on  simulation  results  from 
this  study  will  be  transitioned  through  the  “Sustainability  Logistics  Basing 
(SLB-STO-D)  Demonstrations  and  through  the  Product  Manager  Force 
Sustainment  System-Program  Executive  Office  “Combat  Support  &  Com¬ 
bat  Service  Support”  as  well  as  through  the  Air  Force  Basic  Expeditionary 
Airfield  Resources  (BEAR)  Program. 
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Acronyms  and  Abbreviations 


Term 

Definition 

ACH 

air  changes  per  hour 

CERL 

Construction  Engineering  Research  Laboratory 

DB 

dry  buib 

DoE 

U.S.  Department  of  Energy 

DoD 

U.S.  Department  of  Defense 

ECU 

Environmentai  Controi  Unit 

EFOB-L 

ERDC-CERL  Forward  Operating  Base  Laboratory 

ELA 

equivaient  ieakage  areas 

EMS 

Energy  Management  System 

ERDC 

Engineer  Research  and  Deveiopment  Center 

ERDC-CERL 

Engineer  Research  and  Deveiopment  Center,  Construction 
Engineering  Research  Laboratory 

HFT 

heat  fiux  transducer 

lECU 

improved  Environmentai  Controi  Unit 

JP8 

Jet  Propuision  8 

MEP 

Mountain  Energy  Partnership 

NREL 

Nationai  Renewabie  Energy  Laboratory 

OECIF 

Operationai  Energy  Capitai  improvement  Fund 

ORNL 

Oak  Ridge  Nationai  Laboratory 

PCM 

phase  change  materiai 

PTAC 

Packaged  Terminai  Air  Conditioner 

RMS 

root  mean  square 
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